Preprint typeset in JHBP style - HYPER VERSION 



Unified Dark Matter Scalar Field Models 



Daniele Bertacca"'^"^, Nicola Bartolo"'^, Sabino Matarrese'* '' 

° Dipartimento di Fisica Galileo Galilei Universitd di Padova, via F. Marzolo, 8 1-35131 
Padova, Italy 

^ INFN Sezione di Padova, via F. Marzolo, 8 1-35131 Padova, Italy 

Institute of Cosmology & Gravitation, University of Portsmouth, Dennis Sciama 
Building, Portsmouth, POl 3FX, United Kingdom 

E-mails: daniele . bertaccaOpd . inf n . it , daniele . bertaccaOport . ac . uk, 
nicola. bartoloOpd. inf n. it, sabino .matarreseOpd. inf n. it 

Abstract: In this work we analyze and review cosmological models in which the dynamics 
of a single scalar field accounts for a unified description of the Dark Matter and Dark En- 
ergy sectors, dubbed Unified Dark Matter (UDM) models. In this framework, we consider 
the general Lagrangian of k- essence, which allows to find solutions around which the scalar 
field describes the desired mixture of Dark Matter and Dark Energy. 

We also discuss static and spherically symmetric solutions of Einstein's equations for a 
scalar field with non-canonical kinetic term, in connection with galactic halo rotation 
curves. 
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1. Introduction 

In the last few decades a standard cosmological "Big Bang" model has emerged, based on 
Einstein's theory of gravity, General Relativity. Indeed, observations tell us that - by and 
large - the Universe looks the same in all directions, and it is assumed to be homogeneous 
on the basis of the "Cosmological Principle", i.e. a cosmological version of the Copernican 
principle. The request for the Universe to be homogeneous and isotropic translates, in 
the language of space-time, in a Robertson- Walker metric. Assuming the latter, Einstein 
equations simplify, becoming the Friedmann equations, and in general the solutions of these 
equations are called Priedmann-Lemaitre-Robertson- Walker (FLRW) models. The cosmo- 
logical inhomogeneities we observe on the largest scales as tiny anisotropies of the Cosmic 
Microwave Background (CMB) are then well explained by small relativistic perturbations 
of these FLRW "background" models, while on smaller scales the inhomogeneities are larger 
and call for non-linear dynamics, but relativistic effects arc negligible and Newtonian dy- 
namics is sufficient to explain the formation of the structures we see, i.e. galaxies, groups 
and clusters forming the observed "cosmic web" . In this context, last decade's observations 
of large scale structure, search for la supernovae (SNIa) [1, 2, 3, 4] and measurements of 
the CMB anisotropies [5, 6] suggest that two dark components govern the dynamics of 
the Universe. They are the dark matter (DM), thought to be the main responsible for 
structure formation, and an additional dark energy (DE) component that is supposed to 
drive the measured cosmic acceleration [7, 8]. However, the DM particles have not yet 
been detected in the lab, although there are hints for their existence from cosmic rays 
experiments [9, 10, 11], and there is no theoretical justification for the tiny cosmological 
constant [12] (or more general DE component [7, 8]) implied by observations (see also [13]). 
Therefore, over the last decade, the search for extended theories of gravity has flourished as 
a possible alternative to DE [7, 8]. At the same time, in the context of General Relativity, 
it is very interesting to study the possibility of an interaction between Dark Matter and 
Dark Energy without violating current observational constraints [7, 8, 14, 15, 16, 17, 18, 19] 
(see also [20]). This possibility could alleviate the so called "coincidence problem", namely, 
why are the energy densities of the two dark components of the same order of magnitude 
today. Another more radical explanation of the observed cosmic acceleration and structure 
formation is to assume the existence of a single dark component: Unified Dark Matter 
(UDM) models, see e.g. [21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 
39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50] (see also [51, 52, 53, 54, 55] on how to unify DM, 
DE, and inflation, [56] on unification of DM and DE in the framework of super symmetry, 
[57, 58, 59, 60] on unification of DM and DE from the solution of the strong CP-problem, 
[61, 62] on unification of DM and DE in connection with chaotic scalar field solutions in 
Friedmann-Robertson- Walker cosmologies, [63, 64, 65] on how to unify dark energy and 
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dark matter through a complex scalar field, and [66, 67, 68] on a study of a scalar field, 
"Cosmos Dark Matter", that induces a time-dependent cosmological constant). 

In comparison with the standard DM + DE models (e.g. even the simplest model, 
with DM and a cosmological constant), these models have the advantage that we can 
describe the dynamics of the Universe with a single scalar field which triggers both the 
accelerated expansion at late times and the LSS formation at earlier times. Specifically, for 
these models, we can use Lagrangians with a non-canonical kinetic term, namely a term 
which is an arbitrary function of the square of the time derivative of the scalar field, in the 
homogeneous and isotropic background. 

Originally this method was proposed to have inflation driven by kinetic energy, called 
/c-inflation [69, 70], to explain early Universe's inflation at high energies. Then this scenario 
was applied to DE [71, 72, 73]. In particular, the analysis was extended to a more general 
Lagrangian [74, 75] and this scenario was called fc-essence (see also [71, 76, 77, 78, 79, 80, 
81, 73, 82, 83, 84, 85, 86]). 

For UDM models, several adiabatic or, equivalently, purely kinetic models have been 
investigated in the literature. For example, the generalised Chaplygin gas [22, 23, 24] 
(see also [87, 25, 27, 88, 89, 90, 29, 28, 91, 92, 93, 94]), the Scherrer [30] and generalised 
Scherrer solutions [33], the single dark perfect fluid with "afflne" 2-parameter barotropic 
equation of state (see [37, 39] and the corresponding scalar field models [36]) and the 
homogeneous scalar field deduced from the galactic halo space-time [95, 34]. In general, in 
order for UDM models to have a background evolution that fits observations and a very 
small speed of sound, a severe fine-tuning of their parameters is necessary (see for example 
[39, 29, 25, 28, 27, 30, 31, 96]). Finally, one could also easily reinterpret UDM models 
based on a scalar field Lagrangian in terms of generally non-adiabatic fluids [97, 98] (see 
also [33, 38]). For these models the effective speed of sound, which remains defined in 
the context of linear perturbation theory, is not the same as the adiabatic speed of sound 
(see [99], [70] and [100]). In [38] a reconstruction technique is devised for the Lagrangian, 
which allows to find models where the effective speed of sound is small enough, such that 
the fc-csscnce scalar field can cluster (see also [41, 46, 48, 49, 50]). 

One of the main issues of these UDM models is whether the single dark fluid is able 
to cluster and produce the cosmic structures we observe in the Universe today. In fact, a 
general feature of UDM models is the appearance of an effective sound speed, which may 
become significantly different from zero during the evolution of the Universe. In general, 
this corresponds to the appearance of a Jeans length (or sound horizon) below which the 
dark fluid does not cluster. Thus, the viability of UDM models strictly depends on the value 
of this effective sound speed [99, 70, 100], which has to be small enough to allow structure 
formation [27, 31, 32] and to reproduce the observed pattern of the CMB temperature 
anisotropies [25, 32]. 

In general, in order for UDM models to have a very small speed of sound and a 
background evolution that fits the observations, a severe fine tuning of their parameters is 
necessary. In order to avoid this fine tuning, alternative models with similar goals have been 
analyzed in the literature. Ref. [44] studied in detail the functional form of the Jeans scale 
in adiabatic UDM perturbations and introduced a class of models with a fast transition 
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between an early Einstein-de Sitter cold DM-like era and a later ACDM-likc phase. If the 
transition is fast enough, these models may exhibit satisfactory structure formation and 
CMB fluctuations, thus presenting a small Jeans length even in the case of a non-negligible 
sound speed. Ref. [45] explored unification of DM and DE in a theory containing a scalar 
field of non-Lagrangian type, obtained by direct insertion of a kinetic term into the energy- 
momentum tensor. Finally, Ref. [47] introduced a class of field theories where comprises 
two scalar fields, one of which is a Lagrange multiplier enforcing a constraint between the 
others field value and derivative in order to have the sound speed is always identically zero 
on all backgrounds. 

This work is organized as follows. In Section 2, considering the general Lagrangian of 
k-essence models, wc layout the basic equations. In Section 3 we present an analytical study 
of the Integrated Sachs- Wolfe (ISW) effect within the framework of UDM. Computing the 
temperature power spectrum of the Cosmic Microwave Background anisotropies one is able 
to isolate those contributions that can potentially lead to strong deviations from the usual 
ISW effect occurring in a ACDM Universe. This helps to highlight the crucial role played 
by the sound speed in the unified dark matter models. Our treatment is completely general 
in that all the results depend only on the speed of sound of the dark component and thus it 
can be applied to a variety of unified models, including those which are not described by a 
scalar field but relics on a single dark fluid; see also [32]. In Section 4 wc study and classify 
UDM models defined by the purely kinetic model. We show that these models have only 
one late-time attractor with equation of state equal to minus one (cosmological constant) . 
Studying all possible solutions near the attractor which describes a unifled dark matter 
fluid; see also [33]. Subsequently, noting that purely kinetic models can be described as 
adiabatic single fluid, for these Lagrangians it is natural to give a graphical description on 
pressure - energy density plane, (see also [44]). In Section 5, we present the simplest case 
of a scalar field with canonical kinetic term which unavoidably leads to an effective sound 
speed equal to the speed of light. In Section 6, making the stronger assumption that the 
scalar field Lagrangian is exactly constant along solutions of the equation of motion, we find 
a general class of k- essence models whose classical trajectories directly describe a unified 
Dark Matter/Dark Energy (cosmological constant) fluid. In particular we consider more 
general models allow for the possibility that the speed of sound is small during Einstein- 
de Sitter CDM-like era. In Section 7, we investigate the class of UDM models studied In 
Ref. [38] , which designed a reconstruction technique of the Lagrangian, allowing one to find 
models where the effective speed of sound is small enough, and the fe-essence scalar field 
can cluster (see also [41, 46, 48, 50]). In particular, the authors of Ref. [38] require that the 
Lagrangian of the scalar field is constant along classical trajectories on cosmological scales, 
in order to obtain a background identical to the background of the ACDM model. In Section 
8, we develop and generalize the approach studied in Ref. [38]. Specifically, we focus on 
scalar-field Lagrangians with non-canonical kinetic term to obtain UDM models that can 
mimic a fluid of dark matter and quintessence-like dark energy, with the aim of studying 
models where the background does not necessarily mimic the ACDM background, see also 
[49] . In Section 9, we investigate the static and spherically symmetric solutions of Einstein's 
equations for a scalar field with non-canonical kinetic term, assumed to provide both the 



-4- 



dark matter and dark energy components of the Universe; see also [34]. We show that there 
exist suitable scalar field Lagrangians that allow to describe the cosmological background 
evolution and the static solutions with a single dark fluid. In Section 10, we draw our main 
conclusions. Finally, in the Appendix A, for completeness we provide the spherical collapse 
top-hat solution for UDM models based on purely kinetic scalar eld Lagrangians, which 
allow us to connect the cosmological solutions to the static congurations. 



2. Unified Dark Matter Scalar field models 



We start recalling the main equations which are useful for the description of most the UDM 
models within the framework of k-essence. 
Consider the action 



S = Sg + S^ = J d'^xy^ 



R 



where 



(2.1) 



(2.2) 



where the symbol V denotes covariant differentiation. We adopt SttG = (? = 1 units and 
the (—,+,+,+) signature for the metric (Greek indices run over spacetime dimensions, 
while Latin indices label spatial coordinates). 

The stress-energy tensor of the scalar field ^p has the following form: 



2 5S^ _ dC{^, X) 



and its equation of motion reads 



dX 



d{d^ip)\ dip 



dc 



(2.3) 



(2.4) 



If X is time-like then S^p describes a perfect fluid T'^v = {p + p)Ufj^Uy +pgfj,u, where the 
pressure is 

jC=p{ip,X), (2.5) 



and the energy density is 



The four-velocity has the following form. 



dx 



^/2X 



(2.6) 



(2.7) 



Assume a flat, homogeneous Priedmann-Lemaitre-Robertson- Walker (FLRW) background 
metric, i.e. 

ds"^ = -dt^ + a{tfdijdx'dx^ = a{'qf{-drf' + 6ijdx'dx^) , (2.8) 
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where a{t) is the scale factor, 5ij denotes the unit tensor and r] is the conformal time. 

Assuming that the energy density of the radiation is neghgible at the times of interest, 
and disregarding also the small baryonic component^, the background evolution of the 
Universe is completely characterised by the following equations: 

= a'^H^ = ^a'^p , (2.9) 

and 

n' -H'^ = a'^H = --a'^{p + p) , (2.10) 

where H = a'/a and H = a/ a. The dot denotes differentiation with respect to (wrt) the 
cosmic time t whereas a prime denotes differentiation wrt the conformal time rj. 

In the background we have that X = = (p'^/{2a^), therefore the equation of 

motion Eq. (2.4) for the homogeneous mode ip{t) becomes 

An important quantity is the Equation of State (EoS) parameter w = p/p, which in our 
case reads 

" 2X{dp/dx) -p ■ ^^'-^^^ 

We mainly focus on the other relevant physical quantity, the speed of sound, which enters 
in governing the evolution of the scalar field perturbations. Consider small inhomogeneities 
of the scalar field, i.e. 

ip{t,x) =^o{t) + Sip{t,x) , (2.13) 
and write the perturbed FLRW metric in the longitudinal gauge as 

ds'^ = -(1 + 2$)dt2 + (1 _ 2^)a{tfSijdx'dx^ , (2.14) 

being 6T- = for i / j [101]. The linearised (0 — 0) and {0 — i) Einstein equations are (see 
Ref. [70] and Ref. [100]) 



xrr'po . dp dp p + p 

d(l) dX cl iWaJ v'o 

and 



+ m{p + p)^-^, (2.15) 



5T-^ = -{p + P){'-^) , (2.16) 

where one defines a "speed of sound" relative to the pressure and energy density fluctu- 
ation of the kinetic term [70] as follows: 



2 _ dp/dX dp/dX 



' dp/dx {dp/dx) + 2X{d'^p/dx^) ' ^^'"^^^ 



^Indeed the density of baryons relative is about 4.5% today and 16.9% prior to Dark Energy domination 
in the standard cosmological model [5, 6]. 
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Prom the above hnearized Einstein's equations one obtains [70, 100] 

2 cin 



V^^ = ^^%^f^^ + H , (2.18) 



and 



2$V la'^ip + p) f^5ip 



a 



nj-2 w K^'fo+V- 

Eqs. (2.18) and (2.19) are sufficient to determine the gravitational potential <5 and the 
perturbation of the scalar field. It is useful to write explicitly the perturbed scalar field as 
a function of the gravitational potential 

^ = 2^^, (2.20) 

Defining two new variables 

where z = a'^{p + p)^^'^ /{cgTi), we can recast (2.18) and (2.19) in terms of u and v [100]: 

CsV^u = z{^y , Csv = e(^y (2.22) 

where 9 = l/{csz) = (1 +p/p)~^/^/(\/3a). Starting from (2.22) we arrive at the following 
second order differential equations for u [100]: 

9" 

u" - c^^V^u - — n = . (2.23) 

Unfortunately, we do not know the exact solution for a generic Lagrangian. However, we 
can consider the asymptotic solutions, i.e. the long-wavelength and the short-wavelength 
perturbations, depending whether Cg/c^ <C \9"/9\ or Cgfe^ 2> \9" /9\, respectively. 
Starting from Eq. (2.23), let us define the squared Jeans wave number [32]: 



k 



2 



c'^9 



(2.24) 



Its reciprocal defines the squared Jeans length: Aj 2± a^/kj. 

There are two regimes of evolution. If k'^ S> fcj and the speed of sound is slowly 
varying, then the solution of Eq. (2.23) is 

C 



u ~ 



: exp (±ik J Csdr]^ , (2.25) 



where C is an appropriate integration constant^. On these scales, smaller than the Jeans 
length, the gravitational potential oscillates and decays in time, with observable effects on 
both the CMB and the matter power spectra [32]. 

^This solution is exact if the speed of sound satisfies the equation 2c" Cs — 3 {c'^)^ = 0, which implies 

4 



where ci and C2 are generic constants. A particular case is when ci = 0, for which the speed of sound is 
constant. 
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For large scale perturbations, when k'^ kj, Eq. (2.23) can be rewritten as u" /u 2± 
9" /9, with general solution 



In this large scale limit the evolution of the gravitational potential $ depends only on the 
background evolution, encoded in 9, i.e. it is the same for all k modes. The first term ki9 
is the usual decaying mode, which we are going to neglect in the following, while K2 is 
related to the power spectrum, see e.g. [100]. 

A general feature of UDM models is the possible appearance of an effective sound 
speed, which may become significantly different from zero during the Universe evolution, 
then corresponding in general to the appearance of a Jeans length (i.e. a sound horizon) 
below which the dark fluid does not cluster (e.g. see [99, 32, 39]). Moreover, the presence 
of a non-negligible speed of sound can modify the evolution of the gravitational potential, 
producing a strong Integrated Sachs Wolfe (ISW) effect [32]. Therefore, in UDM models 
it is crucial to study the evolution of the effective speed of sound and that of the Jeans 
length. In other words, one would conclude that any UDM model should satisfy the 
condition that fcj S> k'^ for all scales of cosmological interest, in turn giving an evolution 
for the gravitational potential $ as in Eq. (2.26): 



where = $k (0) Tm (k), $k (0) is the primordial gravitational potential at large scales, 
set during inflation, and (k) is the matter transfer function, sec e.g. [102]. 

Therefore the speed of sound plays a major role in the evolution of the scalar field 
perturbations and in the growth of the over-densities. If Cg is significantly different from 
zero it can alter the evolution of density of linear and non-linear perturbations [99] . When 
cs becomes large at late times, this leads to strong deviations from the usual ISW effect of 
ACDM models [32]. 

In the next section we will perform an analytical study of the Integrated Sachs- Wolfe 
(ISW) effect within the framework of Unified Dark Matter models based on a scalar field 
which aim at a unified description of dark energy and dark matter. Computing the angular 
power spectrum of the Cosmic Microwave Background temperature anisotropics we are able 
to isolate those contributions that can potentially lead to strong deviations from the usual 
ISW effect occurring in a ACDM universe. This helps to highlight the crucial role played 
by the sound speed in the unified dark matter models. 

3. Analytical approach to the ISW effect 

In this Section we focus on the contribution to the large-scale CMB anisotropics which is due 
to the evolution in time of the gravitational potential from the epoch of last scattering up 
to now, the so called late Integrated Sachs- Wolfe (ISW) effect [103]. Through an analytical 
approach we point out the crucial role of the speed of sound in the unified dark matter 
models in determining strong deviations from the usual standard ISW occurring in the 




(2.26) 




(2.27) 
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ACDM models. Our treatment is completely general in that all the results depend only on 
the speed of sound of the dark component and thus it can be applied to a variety of models, 
including those which are not described by a scalar field but relies on a single perfect dark 
fluid. In the case of ACDM models the ISW is dictated by the background evolution, which 
causes the late time decay of the gravitational potential when the cosmological constant 
starts to dominate [104]. In the case of the unified models there is an important aspect to 
consider: from the last scattering to the present epoch, the energy density of the Universe 
is dominated by a single dark fiuid, and therefore the gravitational potential evolution is 
determined by the background and the perturbation evolution of just such a fluid. As a 
result the general trend is the appearance of a sound speed signiflcantly different from zero 
at late times corresponding to the appearance of a Jeans length (or a sound horizon) under 
which the dark fluid does not cluster any more, causing a strong evolution in time of the 
gravitational potential (which starts to oscillate and decay) and thus a strong ISW effect. 
Our results show explicitly that the CMB temperature power spectrum Q for the ISW 
effect contains some terms depending on the speed of sound which give a high contribution 
along a wide range of multipoles i. As the most straightforward way to avoid these critical 
terms one can require the sound speed to be always very close to zero. Moreover we find 
that such strong imprints from the ISW effect come primarily from the evolution of the 
dark component perturbations, rather than from the background expansion history. 
The ISW contribution to the CMB power spectrum is given by 

2'+l„ISW_ 1 /"~<!*:,3|eP"('fo.*)|' 



where Qj^^ is the fractional temperature perturbation due to ISW effect 



ep^(r?o,fe) 

21 + 1 



/ ^'(fi,k)ji[k{m-fl)]dfl, (3.2) 



with r/o and ry* the present and the last scattering conformal times respectively and ji are the 
spherical Bessel functions. Let us now evaluate analytically the power spectrum (3.1). As a 
first step, following the same procedure of Ref. [104], we notice that, when the acceleration 

of the Universe begins to be important, the expansion time scale = i^iw = —1/2) sets 
a critical wavelength corresponding to A;r/i/2 = 1- It is easy to see that if we consider the 
ACDM model then 7/1/2 = "^A i-e. when a^/ao = {^Io/^a)^^^ [10^]. Thus at this critical 
point we can break the integral (3.1) in two parts [104] 



where 



and 



47r ' 27r2 

|2 



Ie,(fc.v.<l) = y^ (3.4) 



/e,(^.V^>l)./^^^^^Jfe --,,--r (3.5) 
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As explained in Ref. [104] the ISW integrals (3.2) takes on different forms in these two 
regimes 

isw(??0,fc) ^ f 2A$fc ji[k{r]o - r?i/2)] A:r/i/2 <Cl . . 

2/ + 1 \ 2$;(%)/i/fc » 1 ^ ' ' 

where A<l>fc is the change in the potential from the matter-dominated (for example at 
recombination) to the present epoch rjo and r]k — i]o — (I + 1/2) /k is the conformal time 
when a given k-mode contributes maximally to the angle that this scale subtends on the 
sky, obtained at the peak of the Bessel function j^. The first limit in Eq. (3.6) is obtained 
by approximating the Bcsscl function as a constant evaluated at the critical epoch 771/2- 
Since it comes from perturbations of wavelengths longer than the distance a photon can 
travel during the the time ^ kick (2A<I>fc) to the photons is the main result, and it will 
corresponds to very low multipolcs, since ryi/2 is very close to the present epoch rjo. It thus 
appears similar to a Sachs- Wolfe effect (or also to the early ISW contribution). The second 
limit in Eq. (3.6) is achieved by considering the strong oscillations of the Bessel functions 
in this regime, and thus evaluating the time derivative of the potentials out of the integral 
at the peak of the Bessel function, leaving the integral [104] 

r-r,^ \/^r[(/ + l)/2] 

With this procedure, replacing (3.6a) in (3.4) and (3.6b) in (3.5) we can obtain the ISW 
contribution to the CMB anisotropics power spectrum (3.1). 

Now we have to calculate, through Eqs. (2.25)-(2.26) and (2.21), the value of ^{k, 77) for 
kili/2 ^ 1 and kr)i/2 3> 1. As we will see that main differences (and the main difficulties) 
of the unified dark matter models with respect to the ACDM case will appear from the 
second regime of Eq. (3.6). 

3.1 Derivation of Iq^ for modes A;77i/2 < 1 

In the UDM models when fc?7i/2 ^ 1 then c^k^ <^ l^'V^I always satisfied. This is due to 
the fact that before the dark fluid starts to behave dominantly as a cosmological constant, 
for rj < r?i/2) its sound speed generically is very close to zero in order to guarantee enough 
structure formation, and moreover the limit krji/2 ^ 1 involves very large scales (since 771/2 
is very close to the present epoch) . For the standard ACDM model the condition is clearly 
satisfied. In this situation we can use the relation (2.26) and can be expressed as in Eq. 
(2.27). The integral in Eq. (2.27) may be written as follows 



/ a^{f})df} = Ir+ I a^{fi)df]. 



(3.8) 



where Ir = a^{fi)dfi and rjR is the conformal time at recombination. When rji < ij < r)R 
the UDM models behave as dark matter In this temporal range the Universe is dom- 
inated by a mixture of "matter" and radiation and Ir = rj^.a^q [{(,r/^) + Cij + (^'^^Z^)]) 

^In fact the Scherrer [30] and generalized Scherrer solutions [33] in the very early Universe, much before 
the equality epoch, have Cs =^ and w > 0. However at these times the dark fluid contribution is sub- 
dominant with respect to the radiation energy density and thus there is no substantial effect on the following 
equations. 
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eg is the value of the scalar factor at matter-radiation equality, = 'q/rj^ and 
^ /24)~^/^ = i^eqlW^ — !)• With these definitions it is easy to see that 



where a, 
V 

o-R = O'eqi^R + 2^r). Notice that Eq. (2.27) is obtained in the case of adiabatic per- 
turbations. Since we are dealing with UDM models based on a scalar field, there will 
always be an intrinsic non-adiabatic pressure (or entropic) perturbation. However for the 

very long wavelengths, ^ 1 under consideration here such an intrinsic perturbation 

turns out to be negligible [70]. For adiabatic perturbations ^k{VR) — (9/10)<I>/fc(0) [101] 
and accounting for the primordial power spectrum, fc^|$/;(0)|^ = Bk"'~^ , where n is the 
scalar spectral index, we get from Eq. (3.6a) 

"V'?l/2 ^ 



^,(fc^i/2<l)~4(2/ + l)S 



Jo 



-k''-'jf[kivo-m/2)] 



1 

To 



a?{rio 



a {fj)dfj 



(3.9) 



where we have neglected //j since it gives a negligible contribution. 

A first comment is in order here. There is a vast class of UDM models that are able to 
reproduce exactly the same background expansion history of the Universe as the ACDM 
model (at least from the recombination epoch on wards) . For such cases it is clear that the 
low i contribution (3.9) to the ISW effect will be the same that is predicted by the ACDM 
model. This is easily explained considering that for such long wavelength perturbations 
the sound speed in fact plays no role. 

3.2 Derivation of Iq^ for modes krji/2 > 1 

As we have already mentioned in the previous section, in general a viable UDM must have 
a sound speed very close to zero for 77 < 77^/2 in order to behave as dark matter also at the 
perturbed level to form the structures we see today, and thus the gravitational potential 
will start to change in time for r] > 771/2- Therefore for the modes krii/2 > 1, in order to 
evaluate Eq. (3.6b) into Eq. (3.5) we can impose that 77^ > 771/2 which, from the definition 
of 77jk 2± 770 — (Z + l/2)/k, moves the lower limit of Eq. (3.5) to (Z + l/2)/{r)o — 771/2). 
Moreover we have that 771/2 ~ ?7o- We can use this property to estimate any observable 
at the value of rjk- Defining x = v/'ni/2 ) and k = A:??i/2, we have = a{r)k) = a(Xfe) = 
ao + {da/dx)\^JXk = 1 - 'ni/2'Ho{l + l/2)/«;, taking oq = 1, and 



dx 



(Xfe) = m/2^'{r]k) 



dx 



XO 



dx^ 



XO 



1 + 1/2 



K 



(3.10) 



where 6xk = Xk — Xo = iVk — %)/^i/2 = ~(^ + ^/2)/n ■ Notice that the expansion (3.10) 
is fully justified, since as already mentioned above, the minimum value of k in Eq. (3.5) 
moves to (Z + l/2)/(7/o/??i/2 — !)> rnaking Sxk much less than 1. Therefore we can write 



l&i isvy(?7o,fc)r 
(21 + 1)2 



= 4 



KiVk)Ii 
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dx^ 



d^k 



dx 



iXk 



d<^i. 



dx 



ixo 



ixo 



1 + 1/2 



— 2r\xo, 



dx' 



1 + 1/2 



(3.11) 
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In this case, during ?/i/2 < ^? < there wiU be perturbation modes whose wavelength stays 
bigger than the Jeans length or smaller than it, i.e. we have to consider both possibilities 
c^k^ <C \0"/0\ and c^k'^ » l^'V^I- general the sound speed can vary with time, and in 
particular it might become significantly different from zero at late times. However, just as 
a first approximation, we exclude the intermediate situation because usually 7/1/2 is very 
close to 770 (see also Ref. [32]). 

3.2.1 Perturbation modes on scales bigger than the Jeans length 

We can see that for n ~ 1 and for Z S> 1 the contribution to the angular power spectrum 
from the modes under consideration is 

^l±}lrisw u, 1 -, /e,(fcr?i/2>l) 1 

In other words we find a similar slope as found for the ACDM model in Refs. [104, 105]. 
Recalling the results of the previous section, this means that in UDM models the con- 
tribution to the ISW effect from those perturbations that are outside the Jeans length 
is very similar to the one produced in a ACDM model. The main difference on these 
scales will be present if the background evolution is different from the one in the ACDM 
model, but for the models where the background evolution is the same, as those proposed 
in Refs. [30, 33, 106, 107, 36] no difference can be observed. 

3.2.2 Perturbation modes on scales smaller than the Jeans length 

When c^k'^ ^ l^'V^i o^*^ must use the solution (2.25) and through the relation (2.21a) the 
gravitational potential is given by 

Mv) = \ [(p + />)/c.]'/' (r?)Cfc(r?i/2)cos ^A:^' c,(7?)d^^ . (3.12) 

In Eq. (3.12) 6^(771/2) = ^*fe(0)Ci/2 is a constant of integration where 



W2 = 2 — ; — , , ,1/2/ ^ ' ^^-^^i 



and it is obtained under the approximation that for 77 < 771/2 one can use the longwavelength 
solution (2.27), since for these epochs the sound speed must be very close to zero. Notice 
that Eq. (3.12) shows clearly that the gravitational potential is oscillating and decaying in 
time. Defining C = C^^^ [{p + p)/ Cg] {ijo) / 4, we take the time derivative of the gravitational 
potential appearing in Eq. (3.6b) by employing the expansion of Eq.(3.11). We thus find 
that, for Cg ~ 1, Eq. (3.5) yields the potentially most dangerous term 



l~ L xo-i 



(3.14) 



with Dq = J^" Cs{x)dx- Such a term makes the angular power spectrum l{l + l)Ci to scale 
as until Z « 25. This angular scale is obtained by considering the peak of the Bessel 
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functions in correspondence of tlie cut-off scale k^q, I ~ keq{f]o — 111/2) ■ fact, for smaller 
scales, l{l + l)Ci will decrease as l/£. This is due to a natural cut-off in the various integrals 
which is introduced for those modes that enter the horizon during the radiation dominated 
epoch, due to the Meszaros effect that the matter fluctuations will suffer until the full 
matter domination epoch. Such a cut-off will show up in the gravitational potential and 
in the various integrals of Eq. (3.14) as a (keq/k)^ factor, where kgq is the wavenumber of 
the Hubble radius at the equality epoch. 

3.3 Discussion of some examples 

Most UDM models have several properties in common. It is easy to see that in Eq. (3.8) 
Ir is negligible because of the low value of Ugq- Moreover in the varioTis models usually 
we have that strong differences with respect to the ISW effect in the ACDM case can be 
produced from perturbations on those scales that are inside the Jeans length as the photons 
pass through them. For these scales the perturbations of the UDM fluid play the main 
role. On larger scales instead we find that they play no role and ISW signatures different 
from the ACDM case can come only from the different background expansion histories. We 
have found that when P > A;} = c'^ \e"/6\ (see (2.23)) one must take care of the term in 
Eq. (3.14). Indeed this term grows faster than the other integrals contained in (3.14) when 
I increases up to ^ ~ 25. It is responsible for a strong ISW effect and hence, in the CMB 
power spectrum /(/ -|- 1)Ci/{2tt), it will cause a decrease in the peak to plateau ratio (once 
the CMB power spectrum is normalized). In order to avoid this effect, a sufficient (but 
not necessary) condition is that the models have satisfy the condition c^k^ < \0"/6\ for the 
scales of interest. The maximum constraint is found in correspondence of the scale at which 
the contribution Eq. (3.14) takes it maximum value, that is A; ~ keq- For example in the 
Generalized Chaplygin Gas model (GCG), i.e when p = — A^/(^+")/p" and = —aw (see 
Section 4), we deduce that \a\ < 10~^ (see Refs. [24] [25] [28]). This is also in agreement 
with the finding of Ref. [27] which performs an analysis on the mass power spectrum and 
gravitational lensing constraints, thus finding a more stringent constraint. 
As far as the generalized Scherrer solution models [33] are concerned, in these models the 
pressure of the UDM fluid is given hy p = gn{X — Xq)"" — A, where gn is a suitable constant 
and n > 1 (see Section 4). The case n = 2 corresponds to unified model proposed by 
Scherrer [30]. In this case we find that imposing the constraint c^fc^ < \9" /6\ for the scales 
of interest we get e = (X - Xq)/Xq < (n - 1) 10"^. 

If we want now to study in greater detail what happens in the GCG model when c^/c^ S> 
l^'Y^I we discover the following things: 

• for 10~^ < a < 5 X 10~^, where we are in the "Intermediate case". Now = —aw 
is very small and the background of the cosmic expansion history of the Universe is 
very similar to the ACDM model. In this situation the pathologies, described before, 
are completely negligible. 

• For 6 X 10-3 < a < 1 a very strong ISW effect is produced; one estimates the same 
order of magnitude for the decrease of the peak to plateau ratio in the anisotropy 
spectrum l{l + l)Q/(27r) (once it is normalized) obtained numerically in Ref. [25] 
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(having assumed that the production of the peaks during the acoustic osciUations at 
recombination is similar to what happens in a ACDM model, since at recombination 
the effects of the sound speed are negligible). 

An important observation arises when considering those UDM models that reproduce the 
same cosmic expansion history of the Universe as the ACDM model. Among these models 
one can impose the condition w = — which, for example, is predicted by UDM models 
with a kinetic term of Born-Infeld type [33, 106, 107, 26]. In this case, computing the 
integral in Eq. (3.14) which gives the main contribution to the ISW effect one can estimate 
that the corresponding decrease of peak to plateau ratio is about one third with respect to 
what we have in the GCG when the value of a is equal to 1. The special case a = 1 is called 
"Chaplygin Gas" (see for example [23]) and it is characterized by a background equation 
of state w which evolves in a different way to the standard ACDM case. Prom these 
considerations we deduce that this specific effect stems only in part from the background 
of the cosmic expansion history of the Universe and that the most relevant contribution to 
the ISW effect is due to the value of the speed of sound c^. 

Let us now make some comments about a particular class of the generalized Chaplygin 
gas models where the sound speed can be larger than the speed of light at late times, i.e. 
when a > 1 (see, for example, [94, 108, 96]). In particular, in [96], the author finds that 
the new constraint a > 350. Indeed, for this range of values, the Jeans wavenumber is 
sufficiently large that the resulting ISW effect is not strong. In this case the Chaplygin 
gas is characterised by a fast transition [44]. However this particular model is ruled out 
because the transition from a pure CDM-like early phase to a post-transition ACDM-like 
late epoch is nearly today {z ~ 0.22). In fact, as discussed in Ref. [44] and in Section 
4.4, the fast transition has to take place sufficiently far in the past. Otherwise, we expect 
that it would be problematic to reproduce the current observations related to the UDM 
parameter w, for instance it would be hard to have a good fit of the CMB and matter 
power spectra. 



4. Purely kinetic Lagrangians 

In this section we focus mainly on Lagrangians C (i.e. the pressure p) that depend only 
on^ X. Defining p{p) = g{X), we have to solve the equation 

p(g{X)) = 2X^^-g{X) (4.1) 
when X is time-like. Then, from Eq. (2.11) we get 

1^ + 2Xp,] ^ + 3 f 2X^] = , (4.2) 
dX dX^ J dN \ dX J ' ^ ' 

where N = In a. We can immediately note that a purely kinetic Lagrangian, through 
Eq. (4.1), (see for example Ref. [33]), can be described as a perfect fluid whose pressure 



''This section is largely based on Ref. [33]. 
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p is uniquely determined by the energy density, since both depend on a single degree of 
freedom, the kinetc term X. In this case = p' / p' corresponds to the usual adiabatic 
sound speed. Obviously if we consider a priori a barotropic or adiabatic equation of state, 
p = p{p) can be described through a purely kinetic k- essence Lagrangian, if the inverse 
function of p = p{p) exists. In Section 4.4, we will use the pressure-density plane to analyze 
the properties that a general barotropic UDM model has to fulfil in order to be viable (see 
also [44]). 

Now we want to make a general study of the attractor solutions in this case. From 
Eq. (2.11) (see Ref. [33]) we obtain the following nodes. 



1) X = X = 0, 2) ^ 



= , (4.3) 

X 



with X a constant. Both cases correspond to w = —1, as one can read from Eq. (2.12). 

In these cases we have either X = or dg/dX = on the node. We know from 
Eq. (4.2) that X can only decrease in time down to its minimum value. This implies that 
w, from Eq. (2.12), will tend to -1 for AT -)• oo. 

At this point wc can study the general solution of the differential equation (4.2). For 
X ^ and dg/dX ^ the solution is [30] 



with k a positive constant. This solution has been also derived, although in a different 
form, in Ref. [109]. As — t- oo, X or dg/dX (or both) must tend to zero, which shows 
that, depending on the specific form of the function g{X), each particular solution will 
converge toward one of the nodes above. From Eq. (4.4), for N ^ oo, the value of X or 
dg/dX (or of both of them) must tend to zero. Then, it is immediate to conclude that 
u; — )■ — 1 is an attractor for N ^ oo and confirms that each of the above solutions will be 
an attractor depending on the specific form of the function g{X). 

In what follows wc will provide some examples of stable node solutions of the equation 
of motion, some of which have been already studied in the literature. The models below 
are classified on the basis of the stable node to which they asymptotically converge. 

4.1 Case 1): Generalized Chaplygin gas 

An example of case 1) is provided by the Generalized Chaplygin (GC) model (see e.g. 
Refs. [22, 23, 24, 25, 28, 29, 110, 27]) whose equation of state has the form 

PGC = -P.[^f , (4.5) 

where now pcc = P and pac = p and p^: and p^ are suitable constants. 

Plugging the equation of state (4.5) into the the continuity equation dpac/dN + 
'^'{PGC + PGc) = 0, we can write pec and pec as function of a. Indeed 



PGC 



-p* 
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(4.6) 
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1-7 



(4.7) 



with u = const. . We note that, when a is smah, we have pcc oc a^"^. In other words, 
this model behaves as DM. Meanwhile, in the late epoch (i.e. a ^ 1), it behaves as a 
cosmological constant. 

Instead, through Eq. (4.1), we can obtain the pressure and the energy density as 
functions of X. Then 



9{X) 



PGC 



-p* 



1/(1-7) 



1-7 



1 

1-7 



1/(1-7) 



1-7 



1-7 



where is a constant. To connect ^ and u we have to use Eq. (4.4). We get 



-p* 



(4.8) 



(4.9) 



(4.10) 



Since = w/^, it is necessary for our scopes to consider the case 7 < 0, so that > 0. 
Note that 7 = — 1 corresponds to the standard "Chaplygin gas" model. Let us obviously 
consider > and > 0. 

Let us conclude this section mentioning two more models that fall into this class of 
solution. The first was proposed in Ref. [Ill], in which g = by/2X — A (with b a suitable 
constant) satisfying the constraint p = —A along the attractor solution Xq = 0. This model, 
however is well-known to imply a diverging speed of sound. The second was proposed 
in Refs. [112, 36, 37, 39] where the single dark perfect fluid with "affine" 2-parameter 
barotropic equation of state p = —A + ap which satisfies the constraint that p = —A along 
the attractor solution Xq = 0. For this model, we have = a, i.e. the speed of sound is 
always a constant. The evolution of p leading to p{a) = A + Pmoo"^'^^'^"^ > where today 
PmO = p(o = 1) — A. When the pressure and the energy density are considered as functions 
of X we have 

g{X) = -A + cX^; p = A + -X^, (4.11) 



where c = pmoct/X^^'^'^^/^'^"^ is the integration constant derived imposing the value of 
the fluid energy density at present and X is X at present time. Prom the matter power 
spectrum constraints [39], it turns out that a < 10~^. 

4.2 Case 2): Scherrer solution 

For the solution of case 1) we want to study the function g around some X = X ^ 0. In 
this case we can approximate 5 as a parabola with ^ lx= 

g = go+92{X-Xf. (4.12) 

with go and g2 suitable constants. This solution, with go < and g2 > 0, coincides with 
the model studied by Scherrer in Ref. [30] (see also Refs. [31, 113]). 
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It is immediate to see that for X — t- X 7^ oo and iV — >■ 00 the value of dX/dN goes to 
zero. Replacing this solution into Eq. (4.4) we obtain 



AglX{X - Xf = ka-^ , 
while the energy density p becomes 

p = -go + ^g2X{X -X) + 3g2{X - Xf 
Now if we impose that today X is close to X so that 

X -X 



e = 



X 



< 1 



then Eq. (4.13) reduces to 



X = X 



1+1^ 
01 



(4.13) 



(4.14) 



(4.15) 



(4.16) 



with ai <C a and with (1/ai) ^ = \l/(2g2)\{k/ X^Y^"^ for e <C 1 . As a consequence, the 
energy density becomes 

.2 Z' « 



P = -50 + 4c/2^ 



ai 



(4.17) 



In order for the density to be positive at late times, we need to impose go < 0. In this case 
the speed of sound (2.17) turns out to be 



{X -X) _1 f a 
{3X -X) ~ 2 \ai 



(4.18) 



We notice also that, for (a/ai)"^ ^ 1 we have ^ 1 for the entire range of validity of 
this solution. Thus, Eq. (4.17) tells us that our k-essence behaves like a fluid with very low 
sound-speed with a background energy density that can be written as 



P = PA + PDM , 



(4.19) 



where pA behaves like a "dark energy" component (p\ = const.) and pdm behaves like a 
"dark matter" component (pdm oc a~^). Note that, from Eq. (4.17), X must be different 
from zero in order for the matter term to be there. (For this particular case the Hubble 
parameter H is a function only of this fluid iJ^ = p/3). 



If the Lagrangian is strictly quadratic in X we can obtain explicit expressions for the 
pressure p and the speed of sound in terms of p, namely 
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(4.20) 
(4.21) 
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Looking at these equations, we observe that in the early Universe {X ^ X i.e. p ^ 
{—go)) the k- essence behaves hke radiation. Therefore, the /c-essence in this case behaves 
Uke a low sound-speed fluid with an energy density which evolves like the sum of a "dark 
matter" (DM) component with p oc and a "dark energy" (DE) component with p = 
const.. The only difference with respect to the standard ACDM model is that in this A;- 
essence model, the dark energy component has Cg <C 1. Starting from the observational 
constraints on /jdm and poE, the value of ai is determined by the fact that the fc-essence 
must begin to behave like dark matter prior to the epoch of matter-radiation equality. 
Therefore, ai < Ugq, where Ogg is the scale factor at the epoch of equal matter and radiation, 
given by = 3 x 10~^ (where we have imposed that the value of the scale factor today 
is ao = 1). At the present time, the component of p corresponding to dark energy in 
equation (4.17) must be roughly twice the component corresponding to dark matter, so 
—go = 852^^(1/01)"^- Substituting ai < Cgg into this equation, we get [30] 



In practice, if we assume that g{X) has a local minimum that can be expanded as a 
quadratic form and when Eq. (4.15) is not satisfied (i.e. for a < ai), we cannot say anything 
about the evolution of X and p. The stronger bound eo < 10~^^ is obtained by Giannakis 
and Hu [31], who considered the small-scale constraint that enough low-mass dark matter 
halos are produced to reionize the Universe. On the other hand the sound speed can be 
made arbitrarily small during the epoch of structure formation by decreasing the value of 
e. One should also consider the usual constraint imposed by primordial nucleosynthesis on 
extra radiation degrees of freedom, which however leads to a weaker constraint. Moreover 
the Scherrer model differs from ACDM in the structure of dark matter halos both because of 
the fact that it behaves as a nearly pressure-less fluid instead of a set of collisioness particles. 
Analytically we will discuss this problem when we will study the static configuration of 
the UDM models, see Section 9 or Ref. [34]. Practically, we will see that when X < 0, 
the energy density of the Scherrer model is negative. Thus, p and p must depend strongly 
on time. In other words, this model will behave necessarily like a fluid and, consequently, 
there is the strong possibility that it can lead to shocks in the non-linear regime [31]. 

4.3 Case 2): Generalized Scherrer solution 

Starting from the condition that we are near the attractor X = X ^ 0, we can generalize 
the definition of g, extending the Scherrer model in the following way 



eo = e(ao = 1) = 



-go 



< 8a^„ < 2 X 10 



1-10 



(4.22) 



p = g = go + gn{x - x) 



(4.23) 



with n > 2 and go and gn suitable constants. 
The density reads 



p = {2n- l)gn{X - X)" + 2Xngn{X - X)"-i - go 
If e" = [{X - X)/XX' < 1 , Eq. (4.4) reduces to 



(4.24) 




(4.25) 
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(where a„_i ^ a) and so p becomes 

p ~ 2nX'^gn - 90 (4.26) 

with (l/a„_i)-^ = [l/(n5f„)](A;/X2"-i)V2for e" < 1. We have therefore obtained the 
important result that this attractor leads exactly to the same terms found in the purely 
kinetic model of Ref. [30], i.e. a cosmological constant and a matter term. One can therefore 
extend the constraint of Ref. [30] to this case, obtaining (eo)"^"'^ = — go / {AnX"^ Qn) < lO^^'^ . 
A stronger constraint would clearly also apply to our model by considering the small-scale 
constraint imposed by the Universe reionization, as in Ref. [31]. If we write the general 
expressions for w and we have 



w = — 



1 + ( ^ ) (X _ X)^ 

90, 



1 - 2nX {X - X)'^-^ - (2n - 1) ( ^] {X - Xf 

\9oJ \goJ 

(4.27) 

c's = ..^^""r^L... (4.28) 



2(n - l)X + (2n - 1){X - X) 
For e ^ 1 one obtains a result similar to that of Ref. [30] , namely 

~ -1 + 2n f ( ' , (4.29) 

VI 50 1/ \an-ij 

On the contrary, when X S> X we obtain 

u; ~ ~ — ^ — (4.31) 
^ 2n-l ^ ' 

In this case one can impose a bound on n so that at early times and/or at high density the 
\i-essence evolves like dark matter. In other words, when n ^ 1, unlike the purely kinetic 
case of Ref. [30], the model is well behaved also at high densities. 

In the section A we study spherical collapse for the generalized Scherrer solution mod- 
els. 

4.4 Studying purely kinetic models in the pressure-density plane. 

In this subsection we report some results of Ref. [44] . Noting that purely kinetic models 
can be described as adiabatic single fluid p = p{p), for these Lagrangians it is natural to 
give a graphical description on the p — p plane, see Fig. 1 . Indeed, this plane gives an idea 
of the cosmological evolution of the dark fluid. Indeed, in an expanding Universe {H > 0) 
Eq. (4.2) implies /? < for a fluid satisfying the null energy condition w > —1 during its 
evolution, hence there exists a one-to-one correspondence between (increasing) time and 
(decreasing) energy density. Finally, in the adiabatic case the effective speed of sound we 
have introduced in Eq. (2.17) can be written as Cg = dp/dp, therefore it has an immediate 
geometric meaning on the p — p plane as the slope of the curve describing the EoS p = p{p)- 
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Figure 1: The UDM p — p plane with the most important areas, (see Ref. [44]). The dashed 
Une represents the p = —p line; the dash-dotted line represents the p = —p/3 line, the boundary 
between the decelerated expansion phase; of the Universe and the accelerated one; the dotted line 
p = —p/10 represents a fictitious boundary, above which the CDM-like behaviour of the UDM fluid 
dominates. The pressure and the energy density are normalised to pA (where pA = A). The ACDM 
model is represented here by the solid horizontal line p/pA = —1, while the line p = represents an 
EdS model, i.e. pure CDM. 



For a fluid, it is quite natural to assume Cg > 0, which then impHes that the function 
p{p) is monotonic, and as such it reaches the p = —p line at some point A.^ Prom the point 
of view of the dynamics this is a crucial fact, because it imphes the existence of an attracting 
fixed point (p = 0) for the conservation equation (4.2) of our UDM fluid, i.e. A plays the 
role of an unavoidable effective cosmological constant. The Universe necessarily evolves 
toward an asymptotic de-Sitter phase, a sort of cosmic no-hair theorem (see [114, 115] and 
refs. therein and [112, 116, 37]).. 

We now summarise, starting from Eqs. (2.9), (2.10) and (4.2) and taking also into 
account the current observational constraints and theoretical understanding, a list of the 
fundamental properties that an adiabatic UDM model has to satisfy in order to be viable. 
We then translate these properties on the p — p plane, see Fig. 1. 

1. We assume the UDM to satisfy the weak energy condition: p ^ 0; therefore, we are 
only interested in the positive half plane. In addition, we assume that the null energy 

® Obviously, we are assuming that during the evolution the EoS allows p to become negative, actually 
violating the strong energy condition, i.e. p < —p/3 at least for some p > 0, otherwise the fluid would never 
be able to produce an accelerated expansion. 
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condition is satisfied: p + p > 0, i.e. our UDM is a standard (non-phantom) fluid. 
Finally, we assume that our UDM models admit a cosmological constant solution 
Lambda at late time, so that an asymptotic equation of state u; = — 1 is built in. 

2. We demand a dust-like behaviour back in the past, at high energies, i.e. a negligible 
pressure p <C p for p A.^ In particular, for an adiabatic fluid we require that at 
recombination l^recl ^ 10~^, see [117, 39, 37, 36]. 

3. Let us consider a Taylor expansion of the UDM EoS p{p) about the present energy 
density po: 

p ~ po + a(p - Po) , (4.32) 

i.e. an "affine" EoS model [39, 37, 36, 112] where a is the adiabatic speed of sound 
at the present time. Clearly, these models would be represented by straight lines in 
Fig. 1, with a the slope. The ACDM model, interpreted as UDM, corresponds to the 
affine model (4.32) with a = (see [112] and [37, 39]) and thus it is represented in 
Fig. 1 by the horizontal line p = —A. From the matter power spectrum constraints 
on afhnc models [39], it turns out that a < 10~^. Note therefore that, from the UDM 
perspective, today we necessarily have w ~ —0.7. 



Few comments are in order. From the points above, one could conclude that any 
adiabatic UDM model, in order to be viable, necessarily has to degenerate into the ACDM 
model, as shown in [27] for the generalised Chaplygin gas and in [39] for the affine adiabatic 
modeF. In other words, one would conclude that any UDM model should satisfy the 
condition <C 1 at all times, so that k!j S> k'^ for all scales of cosmological interest, in 
turn giving an evolution for the gravitational potential $ as in Eq. (2.26). 

On the other hand, let us write down the explicit form of the Jeans wave-number: 



k 



3 p (l + w) 



2 

^ 2(1 + Z)2 c2 



_ _ , 3{ct-wr-2{ci-w) 1 

2^« ' ^dp^ Q{l + w) ^3 



(4.33) 



Clearly, we can obtain a large fcj not only when Cg — > 0, but also when Cg changes rapidly, 
i.e. when the above expression is dominated by the p dc^/dp term. When this term is 
dominating in Eq. (4.33), we may say that the EoS is characterised by a fast transition. 

In the paper [44] the authors investigate observational constraints on UDM models 
with fast transition, introducing and discussing a toy model. In particular, they explore 
which values of the parameters of such a toy model fit the observed CMB and matter power 
spectra. 



®Note that we could have p ~ —A and yet, if p A, the Universe would still be in a matter-like era. 

^From the point of view of the analysis of models in the p — p plane of Fig. 1, the constraints found by 
Sandvik et al [27] on the generalised Chaplygin gas UDM models and by [39] on the affine UDM models 
simply amount to say that the curves representing these models are indistinguishable from the horizontal 
ACDM line. 
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5. UDM Scalar Field with canonical kinetic term 



Starting from the barotropic equation of state p = p{p) we can describe the system either 
through a purely kinetic k- essence Lagrangian, as we already explained in the last section, 
or through a Lagrangian with canonical kinetic term, as in quintessence-like models (see 
Ref. [33]). 

In the second case we have to solve the two differential equations 



X + V{^) = p{^,X) 



(5.1) 
(5.2) 



where X = </3^/2 is time-like. In particular, if we assume that our model describes a unified 
dark matter/dark energy fluid, we can proceed as follows: starting from p = —3H{p + p) = 
—\/3p{p + p) and 2X = {jp + p) = {dip / dp)^ (? we get 



if 



±- 



1 



dp'l^ 



VsJ., {p{p') + p'V/' ' 



(5.3) 



up to an additive constant which can be dropped without any loss of generality. Inverting 
the Eq. (5.3) i.e. writing p = p{ip) we are able to get V{ip) = [p{ip) — p{p{ip))]/2. Now we 
require that the fluid has constant pressure p = —A, i.e. that the Lagrangian of the scalar 
field is constant along the classical trajectory corresponding to perfect fluid behavior. In 
other words one arrives at an exact solution with potential 



+ 1 



(5.4) 



see also Refs. [106, 107]. For large values of 99, V{(p) oc exp(-v/3'/?) (cquivalently, for large 
values of —tp, V{(p) oc cxp{~^/3ip)) and our scalar field behaves just like a pressurelcss dark 
matter fluid. Indeed, this asymptotic form, in the presence of an extra radiation component, 
allows to recover one of the stable nodes obtained in Ref. [118] for quintessence fields with 
exponential potentials, where the scalar field mimics a pressurelcss fluid. Under the latter 
hypothesis we immediately obtain 



V3 



arccosh (p/A)^^"^ 



(5.5) 



which can be inverted to give the scalar field potential of Eq. (5.4) as V{ip) = {p{ip) + A)/2. 
One then obtains 



If = — -\/Asinh f 



(5.6) 



which can be immediately integrated, to give 



C = exp 



{t - Q 



(5.7) 
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for t > t*, with such that ip{t — t- t*) — t- oo. Replacing this solution in the expression for 
the energy density one can easily solve the the Friedmann equation for the scale-factor as 
a function of cosmic time, 



a{t) = ao —TT^ r ' (^-S) 



sinh^/^ 




sinh2/3 





which coincides with the standard expression for a fiat, matter plus Lambda model [119], 
with fioA/^Om = sinh^[^^^(to — t*)], ^OA and ^Om being the cosmological constant and 
matter density parameters, respectively. 

Using standard criteria (e.g. Ref. [8]) it is immediate to verify that the above trajectory 
corresponds to a stable node even in the presence of an extra-fluid (e.g. radiation) with 
equation of state w^uid = Pfluid/Pfluid > 0, where Pfluid and pfluid are the fluid pressure 
and energy density, respectively. Along the above attractor trajectory our scalar field 
behaves precisely like a mixture of pressureless matter and cosmological constant. Using 
the expressions for the energy density and the pressure we immediately find, for the matter 
energy density 

= P — A = Asinh^ ^^'^^ ' ^^'^^ 

The peculiarity of this model is that the matter component appears as a simple consequence 
of having assumed the constancy of the Lagrangian. 

A closely related solution was found by Salopek &; Stewart [120], using the Hamiltonian 
formalism. 

To conclude this section, let us stress that, like any scalar field with canonical kinetic term 
[121, 73], our UDM model predicts = 1, as it is clear from Eq. (2.17), which inhibits the 
growth of matter inhomogcncities. In summary, we have obtained a "quartessence" model 
which behaves exactly like a mixture of dark matter and dark energy along the attractor 
solution, whose matter sector, however is unable to cluster on sub-horizon scales (at least 
as long as linear perturbations are considered). 

6. UDM Scalar Field with non-canonical kinetic term 

We can summarize our findings so far by stating that purely kinetic k- essence cannot 
produce a model which exactly describes a unified fluid of dark matter and cosmological 
constant, while scalar field models with canonical kinetic term, while containing such an 
exact description, unavoidably lead to = 1, in conflict with cosmological structure for- 
mation. In order to find an exact UDM model with acceptable speed of sound we consider 
more general scalar field Lagrangians (see Ref. [33]). 

6.1 Lagrangians of the type jC{ip,X) = g{X) — V{ip) 

Let us consider Lagrangians with non-canonical kinetic term and a potential term, in the 
form 

£{^,X)=g{X)-V{ip) . (6.1) 
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The energy density then reads 



(6.2) 



while the speed of sound keeps the form of Eq. (2.17). The equation of motion for the 
homogeneous mode reads 



dX + ^^dX^ dN + ^ r^dX j - ~dN 



One immediately finds 



One can rewrite the equation of motion Eq. (6.3) in the form 

-- . 



dT 

2X 

dX 



dX f ^ dV 



(6.3) 



(6.4) 



(6.5) 



It is easy to see that this equation admits 2 nodes, namely: 

1) dg/dX\^ = and 

2) X = 0. 

In all cases, for iV — )■ oo, the potential V should tend to a constant, while the kinetic term 
can be written around the attractor in the form 



g{X) = 



X -X 

M4 



n > 2 , 



(6.6) 



where M is a suitable mass-scale and X a constant. The trivial case g{X) = X obviously 
reduces to the one of Section 4. 

Following the same procedure adopted in the previous section we impose the constraint 
p = —A, which yields the general solution pm = 2J^{X). 

This allows to define (p = ip{pm) as a solution of the diflFerential equation 



2J" 



(6.7) 



As found in the case of k- essence, the most interesting behavior corresponds to the 
limit of large n and X = in Eq. (6.6), for which we obtain 



A sinh 



-2 



3A 
8M4 



1/2 



(6.8) 



leading to V{tp) ^ Prnf^n — A, and = l/(2n — 1) ?s 0. The Lagrangian of this model is 
similar to that analyzed in Ref. [95]. 
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6.2 Lagrangians of the type C{^p,X) = f{(p)g{X) 

Let us now consider Lagrangians with a non-canonical kinetic term of the form C((p,X) = 
f{ip)g{X) (see Ref. [33]). 

Imposing the constraint p = —A, one obtains /(v) = 9{^)^ which inserted in the 
equation of motion yields the general solution 

- d\-n.\g\ _ pm 
dX ~ 2A ■ 



X- 



(6.9) 



The latter equation, together with Eq. (6.7) define our general prescription to get UDM 
models describing both DM and cosmological constant-like DE. 

As an example of the general law in Eq. (6.9) let us consider an explicit solution. 
Assuming that the kinetic term is of Born-Infeld type, as in Refs. [26, 122, 106, 107], 

g{X) = -V1-2X/M4 , (6.10) 

with M a suitable mass-scale, which implies p = f{(p)/^/l — 2X/M^, we get 



X{a) 



(6.11) 



2 1 + ka-^ ' 

where k = pm{a*)al/ A and a* is the scale- factor at a generic time t*. In order to obtain 
an expression for 'p{a), we impose that the Universe is dominated by our UDM fluid, i.e. 
= p/3. This gives 

2M2 



3A 



arctan 



(fca-3) 



which, replaced in our initial ansatz p 
Ref. [106, 107]) 



cos 



(6.12) 

A allows to obtain the expression (see also 
^ (6.13) 



If one expands /(<^) around (p = 0, and X/M'^ <C 1 one gets the approximate La- 
grangian 



C 



A 



2M4 



A 



1 + 



3A 
8AF 



(6.14) 



Note that our Lagrangian depends only on the combination (f/M'^, so that one is free to 
reabsorb a change of the mass-scale in the definition of the filed variable. Without any loss 
of generality we can then set M = A^/^, so that the kinetic term takes the canonical form 
in the limit X <C 1. We can then rewrite our Lagrangian as 



^1 - 2X/A 



cos 



(4 



(6.15) 



This model implies that for values of -\/3(^ — tt and 2X/A 1, 



cos 



■3/2 



(6.16) 
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Figure 2: Evolution of the scalar field density parameter vs. redshift (seeRef. [33]). The continuous 
line shows the UDM density parameter; the dashed line is the density parameter of the DM + DE 
components in a standard ACDM model; the dotted line is the radiation density parameter. 



the scalar field mimics a dark matter fluid. In this regime the effective speed of sound is 
= 1 - 2X/A Ri 0, as desired. 

To understand whether our scalar field model gives rise to a cosmologically viable UDM 
solution, we need to check if in a Universe filled with a scalar field with Lagrangian (6.15), 
plus a background fluid of e.g. radiation, the system displays the desired solution where 
the scalar field mimics both the DM and DE components. Notice that the model does not 
contain any free parameter to specify the present content of the Universe. This implies 
that the relative amounts of DM and DE that characterize the present Universe are fully 
determined by the value of (po = ip{to)- In other words, to reproduce the present Universe, 
one has to tune the value of f{(p) in the early Universe. However, a numerical analysis 
shows that, once the initial value of (p is fixed, there is still a large basin of attraction in 
terms of the initial value of dip/dt, which can take any value such that 2X/A <^ 1. 

The results of a numerical integration of our system including scalar field and radiation 
are shown in Figures 2-4. Figure 2 shows the density parameter, Oudm as a function of 
redshift, having chosen the initial value of if so that today the scalar field reproduces the 
observed values JIdm and JIde- Notice that the time evolution of the scalar field energy 
density is practically indistinguishable from that of a standard DM plus Lambda (ACDM) 
model with the same relative abundances today. Figure 3 shows the evolution equation of 
state parameter wubm] once again the behavior of our model is almost identical to that of 
a standard ACDM model for 1 + 2; < 10^. Notice that, since = — u^udm, the effective 
speed of sound of our model is close to zero, as long as matter dominates, as required. In 
Figure 4 we finally show the redshift evolution of the scalar field variables X = 0^/2 and 
ip: one can easily check that the evolution of both quantities is accurately described by the 
analytical solutions above, Eqs. (6.11) and (6.12), respectively (the latter being obviously 
valid only after the epoch of matter-radiation equality). 
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Figure 3: The redshift evolution of the scalar field equation of state parameter wudm (continuous 
line) is compared with that of the sum of the DM + DE components in a standard ACDM model 
(dashed line), see Ref. [33]. 

A 

y 

X 

/V 

T 





2 V3" 



-©■ 

7T 

4 3 2 1 

Figure 4: Redshift evolution of the scalar field of the scalar field variables X = ip^ /2 (top) and ip 
(bottom), see Ref. [33]. 

However in this model, as discussed in Ref. [32], the non-negligible value of the sound 
speed today gives a strong contribution to the ISW effect and produces an incorrect ratio be- 
tween the first peak and the plateau of the CMB anisotropy power-spectrum 1{1+1)Ci/{2'it). 

7. How the Scalar Field of Unified Dark Matter Models Can Cluster 

The authors of [38] proposed a technique for constructing UDM models where the scalar 
field can have a sound speed small enough to allow for structure formation and to avoid a 
strong integrated Sachs- Wolfe effect in the CMB anisotropics which typically plague UDM 
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models*^ (sec also [41, 46]). In particular, they studied a class of UDM models where, at 
all cosmic times, the sound speed is small enough that cosmic structure can form. To do 
so, a possible approach is to consider a scalar field Lagrangian L of the form 

L = pi^, X) = f{<p)g{X) - V{ip) . (7.1) 

Therefore, by introducing the two potentials f{ip) and V{ip), we want to decouple the 
equation of state parameter w and the sound speed Cg. This condition does not occur 
when we consider either Lagrangians with purely kinetic terms or Lagrangians like C = 
g{X) — V{(p) OT C = f{(p)g{X) (see for example [33] and the previous Sections 6.1 and 
6.2). In the following subsections we will describe how to construct UDM models based on 
Eq. (7.1), following the analysis of Ref. [38]. 

7.1 How to construct UDM models 

Let us consider the scalar field Lagrangian of Eq. (7.1). The energy density p, the equation 
of state w and the speed of sound are 



p(x,v) = /(v) 



dX ' 



V{^) , (7.2) 



,y f{<p)g{X)-V{^) 

f{<p) [2X {dg{X)/dX) - g{X)] - V{^) ' ^'"^^ 

c\x) = idg{X)/dX) 

' {dgiX)/dX)+2X{d^g{X)/dX^) ' ^ 

respectively. The equation of motion (2.11) becomes 

dg ,d^g\dX ^dg dlnff^^^dq \ \ dV 



Unlike in models with a Lagrangian with purely kinetic terms, here we have one more 
degree of freedom, the scalar field configuration itself. This allows to impose a new condition 
to the solutions of the equation of motion. In Ref. [33], the scalar field Lagrangian was 
required to be constant along the classical trajectories. Specifically, by requiring that 
£ = — A on cosmological scales, the background is identical to the background of ACDM. 
In general this is always true. In fact, if we consider Eq. (2.11) or, equivalently, the 
continuity equations {dp/dN) = — 3(p + p), and if we impose p = —A, we easily get 

P = PDM(a = 1) a"^ + A = pDM + PA , (7.6) 

where p\ behaves like a cosmological constant "Dark Energy" component (pA = const.) 
and pdm behaves like a "Dark Matter" component (pdm a~^). This result implies that 
we can think the stress-energy tensor of our scalar field as being made of two components: 
one behaving like a pressure-less fluid, and the other having negative pressure. In this way 
the integration constant Pdm(« = 1) can be interpreted as the "dark matter" component 



*This section is largely based on Ref. [38]. 
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today; consequently, flm{0) = Pdm{ci = l)/(3i^^(a. = 1)) and riA(O) = A/(3iI^(a = 1)) are 
the density parameters of "dark matter" and "dark energy" today. 

Let us now describe the procedure that we will use in order to find UDM models 
with a small speed of sound. By imposing the condition jC{X, ip) = —A, we constrain the 
solution of the equation of motion to live on a particular manifold A^a embedded in the 
four dimensional space-time. This enables us to define 99 as a function of X along the 
classical trajectories, i.e. 99 = C^^{X, A)\j^^. Notice that therefore, by using Eq.(7.5) and 
imposing the constraint p = —A, i.e. V{ip) = f{ip)g{X) + A, we can obtain the following 
general solution of the equation of motion on the manifold A^a 

2X^^f{^{X)) =Kua-\ (7.7) 

where u = r2„i(0)/r2A(0) . Here we have constrained the pressure to bep = —A. In Section 
8 we will describe an even more general technique to reconstruct UDM models where the 
pressure is a free function of the scale factor a. 

If we define the function g{X), we immediately know the functional form of with 
respect to X (see Eq. (7.4)). Therefore, if we have a Lagrangian of the type C = f{ip)g{X) 
or £ = g{X) — V{ip), we are unable to decide the evolution of c^{X) along the solutions 
of the equation of motion [33] because, once g{X) is chosen, the constraint C = —A fixes 
immediateley the value of /(<^) or V{ip). On the contrary, in the case of Eq. (7.1), we can 
do it through the function f{ip{X)). In fact, by properly defining the value of f{(p{X)) 
and using Eq.(7.5), we are able to fix the slope of X and, consequently (through g{X)), 
the trend of c^(X) as a function of the scale factor a. 

Finally, we want to emphasize that this approach is only a method to reconstruct the 
explicit form of the Lagrangian (7.1), namely to separate the two variables X and (p into 
the functions g, f and V. 

Let us now give an example where we apply this prescription. In particular, in the 
following subsection, we assume a kinetic term of Born-Infeld type [26, 122, 106, 107]. 
Other examples (where we have the kinetic term g{X) of the Scherrer model [30] or where 
we consider the generalized Scherrer solutions [33] ) are reported in Ref. [38] . 

7.1.1 Lagrangians with Born-Infeld type kinetic term 

Let us consider the following kinetic term 



g{X) = -Vl - 2X/M^ , (7.8) 
with M a suitable mass scale. We get 

2X/M^ /(V9(X)) = A a-^ (7.9) 

and 

cl{X) = 1 - 2A:/M^ . (7.10) 

In the next subsection, we give a Lagrangian where the sound speed can be small. It is 
important to emphasize that the models described here and in the next subsection satisfy 
the weak energy conditions p> and p + p > 0. 
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7.2 UDM models with Born-Infeld type kinetic term and a low speed of sound 

Let us consider for / the following definition 



A 



mx)) = - 



2X/M^ - h 



(1 - 2X/M4)V2 ' 



(7.11) 



where h and jj, are appropriate positive constants. Moreover, we impose that ^ < 1. Thus 
we get 



^, , M^h + jjva-^ 
Xia) = ——- 5- or 



dip 
dN 



h + fih'a 



and, for c^, we obtain the following relation 



3M^ 



A (1 + i^a-3) (1 + /ii/a- 



1 - h 



1 + ijiva 



(7.12) 



(7.13) 



Therefore, with the definition (7.11) and using the freedom in choosing the value of h, we 
can shift the value of for a — t- oo. Specifically, h = 1 — where c^o = Cs{a — t- oo). At 
this point, by considering the case where h = ix (which makes the equation analytically 
integrable), we can immediately obtain the trajectory (/j(a), namely 



ip{a) 



a)=( 



V 3A 



4\ 1/2 



arcsinh(z/^a ^) 



-3\-V2 



(7.14) 



Finally, we obtain 
and 



A(l - /i)V2 cosh [ 


i4hMi) ^_ 




^ sinh 


\( 3A ^V2 1 


|l + /isinh2 


( V^^J 
VAhMi) ^\ 


} 



(7.15) 



V{ip) 



A' 




\( 3A ^^2 1 
[[4hM^) ^\ 


+ 2/i- l} 
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1 + h sinh^ 


'( 3A 
U/iM4 





(7.16) 



This result implies that in the early universe -^3 A/ (4/iM*^) </? <C 1 and 2X/M^ « 1, and 
we obtain 



A/iM^v^^ A^/r^ 1 



oc a 



3/2 



A(2/i- 1) 



h 



\g{X)\ = Vl - 2X/A (X a-=^/2 , 
« /m(2X^^-5(X)] oc a-%7.17) 



In other words, we find, for f{ip) and g{X), a behaviour similar to that we have studied in 
Section 6.2, as also in Ref. [33]. 

When a — >■ oo, we have <^ — >■ oo and 2X/M^ h. Therefore 



f{^)g{X) 



A 
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Figure 5: Sound speed Cs^(a) for different values of Cqo^ = 10~^, 10~^, 10~*, 10"'^ from bottom to 
top (see [41]). 

that is, for a — t- oo, the dark fluid of this UDM model will converge to a Cosmological 
Constant. 

In Ref. [38] the authors analytically show that, once the initial value of is fixed, there 
is still a large basin of attraction in terms of the initial value of dip/dt, which can take any 
value such that 2X/M^ <C 1. Moreover, Ref. [38] investigates the kinematic behavior of 
this UDM fluid during the radiation-dominated epoch. 

We can conclude that, once it is constrained to yield the same background evolution 
as ACDM and we set an appropriate value of Coo, this UDM model provides a sound 
speed small enough that i) the dark fluid can cluster and ii) the Integrated Sachs- Wolfe 
contribution to the CMB anisotropies is compatible with observations. Figure 5 shows an 
example of the dependence of cj. on CL for different values of Cqo • 

In Fig. 6 we present some Fourier components <I>fc(a) of the gravitational potential, 
normalized to unity at early times (see [41]). As we can note from this figure, the possible 
appearance of a sound speed significantly different from zero at late times corresponds to 
the appearance of a Jeans' length under which the dark fluid does not cluster any more, 
causing a strong evolution in time of the gravitational potential. By increasing the sound 
speed, the potential starts to decay earlier in time, oscillating around zero afterwards. 
Moreover, at small scales, if the sound speed is small enough, UDM reproduces ACDM. 
This reflects the dependence of the gravitational potential on the effective Jeans' length 
Aj(a) [32]. 

Finally, in Ref. [41] the authors show for this UDM model the lensing signal in linear 
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Figure 6: Normalized potentials ^k{a)/^kiO) are shown for ACDM (solid) and UDM (dot-dashed) 
(see [41]). The lower panel shows potentials at fc = 0.001 ft- Mpc^^, the medium panel at fc = 
0.01 /iMpc"^ and the upper panel at A: = 0.1 ftMpc^^ UDM curves are for c^^ = 10'^, 10'^, 10"^ 
from top to bottom, respectively. At small Coo, ACDM and UDM curves are indistinguishable. 



theory as produced in ACDM and UDM; as sources, they consider the CMB and back- 
ground galaxies, with different values of the peak and different shapes of their redshift 
distribution. For sound speed lower than Cqo = 10^^, in the window of multipoles / > 10 
(Limber's approximation) and where our ignorance on non-linear effects due to small scales 
dynamics become relevant, the power spectra of the cosmic convergence (or shear) in the 
flat-sky approximation in UDM and ACDM are similar. When the Jeans' length Aj(a) 
increases, the Newtonian potential starts to decay earlier in time (at a fixed scale), or at 
greater scales (at a fixed epoch). This behaviour reflects on weak lensing by suppressing the 
convergence power spectra at high multipoles. They find that, for values of the sound speed 
between Cqo = 10~^ and Cqo = 10^^, UDM models are still comparable with ACDM, while 
for higher values of Coo these models are ruled out because of the inhibition of structure 
formation. Moreover, they find that the dependence of the UDM weak lensing signal on 
the sound speed Cqo increases with decreasing redshift of the sources. They also show the 
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errors for the fiducial ACDM signal for wide-field surveys like EUCLID or Pan-STARRS, 
and they find that one isin principle able to distinguish ACDM from UDM models when 
Coo ^ 10~^. Moreover, in Ref. [46] the authors calculate the 3D shear matrix C'^"i{ki,k2\(.) 
in the flat-sky approximation for a large number of values of Cqo- They see that, whilst the 
agreement with the ACDM model is good for small values of Cqo, when one increases the 
sound speed parameter, the lensing signal appears more suppressed at small scales, and 
the 3D shear matrix shows bumps related to the oscillations of the gravitational potential. 
Moreover, they show that the expected evidence clearly shows that the survey data would 
unquestionably favour UDM models over the standard ACDM model, if its sound speed 
parameter exceed 10~^. 

7.3 Prescription for UDM Models with a generic kinetic term 

We now describe a general prescription to obtain a collection of models that reproduce a 
background similar to ACDM and have a suitable sound speed. Some comments about 
the master equation (7.7) are first necessary. The relation (7.7) enables to determine a 
connection between the scalar factor a and the kinetic term X on the manifold M.\ and 
therefore a mapping between the cosmic time and the manifold M.k. 

Now it is easy to see that the LHS of Eq. (7.7), seen as a single function of X, must 
have at least a vertical asymptote and a zero, and the function must be continuous between 
the two. In particular, when X is near the vertical asymptote the universe approaches the 
cosmological constant regime, whereas when X is close to the zero of the function, the dark 
fluid behaves like dark matter. Therefore, if we define 

«^(^» = 2x(amiax) 

where, for example, 

HX) = , (7.19) 

(where jx is an appropriate positive constant) the value of Xf and Xi are the zero and the 
asymptote mentioned above, namely, when a — t- we have X ^ Xi and when a — t- oo we 
have X — )• Xf. Moreover, if Xf > Xi we have dX/dN > 0, whereas if Xf < Xi we have 
dX/dN < 0. In other words, according to Eq.(7.7), 

X(a)^X,l±^W^. (T.20) 

Let us emphasize that the values of Xi and Xf are very important because they auto- 
matically set the range of values that the sound speed can assume at the various cosmic 
epochs. 

Let us finally make another important comment. One can use this reconstruction of 

the UDM model in the opposite way. In fact, by imposing a cosmological background 
identical to ACDM, the observed CMB power spectrum, and the observed evolution of 
cosmic structures, one can derive the evolution of the sound speed vs. cosmic time. In 
this case, by assuming an appropriate kinetic term g{X) through Eq. (7.4), we can derive 
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X{a) and, consequently, '/^(a) and X(a{{p)) = X{{p). Therefore, by using the relations (7.7) 
and V{ip) = f{ip)g{X) + A, one can determine the functional form of /(</?) and V{ip). 

8. Generalized UDM Models 

In this Section we consider several possible generalizations of the technique introduced in 
Section 7.1, with the aim of studying models where the background does not necessarily 
mimic the ACDM background (sec [38, 49]). 

Let us consider a scalar field Lagrangian C of the form 

CiX, ^) = p{^, X) = f{ip)g{h{ip)X) - V{ip) . (8.1) 

Note that, introducing the three potentials f{ip), h{ip) and V{ip), we follow an approach 
similar to the one studied in Ref. [38] in order to decouple the equation of state parameter 
w and the sound speed Cg. In order to reconstruct these potentials we need three dynamical 
conditions: a) a choice for p{N), b) the continuity equation or, equivalently, the equation 
of motion (2.11), c) a choice for c^{N) (see [49]). 

Let us obtain the Lagrangian through two different simple approaches: 



1) By choosing p{N). Indeed we get 
dp 



dN 



+ 3p = -3p{N) , i.e. p{N) = e 



-3N 



/N 
(^e^'''p{N')dN') + K 



(8.2) 

where K is an integration constant. By imposing the condition J0,{X, cp) = p{N) along 
the classical trajectories, we obtain (p = C~^{X{N),p{N))\ ^ . Thus, starting from 
a generic Lagrangian £, = f{ip)g{h{tp)X) — V{tp) we get 

-dg{h{p^[X,N])X) 



2X{N) 



dx 



{N)f{p{X,N))=p{N)+e- 



3N 



/N 
(e3^'p(iV')diV') + K 

(8.3) 

p{a = 1) — A. The freedom provided by the choice of K 



For example, if p = —A, K 
is particularly relevant. In fact, by setting i^T = 0, we can remove the term p oc a~^. 
Alternatively, when if 7^ 0, we always have a term that behaves like presseure-less 
matter. We thus show that the single fluid of UDM models can mimic not only a 
cosmological constant but also any quintessence fluid. 

Thus, using Eq. (8.3) and following the procedure described in Section 7.1, one gets 
the relations X = Qp{N), and consequently 




/xl/2 



/N 
(e^'''p{N')dN')+Ke- 



-3N 



-1/2 



dN'} .(8.4) 
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Therefore, with the functions Gp{N) and Qp{N), one can write f{X, N) = f{gp{N),N) 
f{Gp{Qp^{^))i Qp^i^)) = /(¥')• Thus, by starting from a Lagrangian whose behav- 
ior is given by p{N), the speed of sound is determined by the appropriate choice of 
g{h{^)X), where h{X,N) = h{gp{N),N) = h{gp{Q-'{ip)), QpH^)) = h{ip). 

2) By choosing the equation of state w{N). Indeed 

p(iV) = poe-=^/"W^W^', (8.5) 

where po is a positive integration constant, and 

p{N) = po^«(Ar)e-3i"'(-(^')+i)'^^' . (8.6) 

Therefore, still by imposing the condition C{X,(p) = p[w{N),N] along the classical 
trajectories, i.e. tp = C~^[X{N),p{w{N), N)] \ , one gets 

^^ dg{h{^[X N])X) ^^^^ ^ ^^^^^^^ ^ i]^-3r (.(iV')-M).iv' . (8.7) 
oX 

Therefore, on the classical trajectory we can impose, by using w{N), a suitable 
function p{N) and thus the function p{N). The master equation Eq. (8.7) generalizes 
Eq. (7.7). Also in this case, by Eq. (8.7) and by following the argument described in 
Section 7.1, one can get the relations X = QwiN), and consequently 



{w{N")+l)dN" 



-1/2 ^ 

dN' }+^o- (8.8) 



Thus, with the functions ^„(iV) and Qw{N), one can write f{X, N) = /((?„, (AQ, = 
fiQwiQw^i^)), Qw^{^)) = /(v)- Then we can find a Lagrangian whose behavior is 
determined by ■w(N) and whose speed of sound is determined by the appropriate 
choice of gih{ip)X), where h{X,N) = h{gp{N),N) = h{gp{Q-\^)),Q-\^)) = 
h{ip). 

Let us conclude that the p{N) constraint on the equation of motion is actually a weaker 
condition than the ■w{N) constraint. The larger freedom that the p(N) constraint provides 
naturally yields an additive term in the energy density that decays like a~^, i.e. like a 
matter term in the homogeneous background. Let us emphasize that this important result 
is a natural consequence of the p{N) constraint and is not imposed a priori (see [38, 49]). 



9. Halos of Unified Dark Matter Scalar Field 

A complete analysis of UDM models should necessarily include the study of static solutions 
of Einstein's field equations. This is complementary to the study of cosmological back- 
ground solutions and would allow to impose further constraints to the Lagrangian of UDM 
models. The authors of Refs. [123] and [95] have studied spherically symmetric and static 
configuration for k-essence models. In particular, they studied models where the rotation 
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velocity becomes flat (at least) at large radii of the halo. In these models the scalar fleld 
pressure is not small compared to the mass-energy density, similarly to what is found in the 
study of general fluids in Refs. [124, 125, 126, 127], and the Einstein's equations of motion 
do not reduce to the equations of Newtonian gravity. Further alternative models have been 
considered, even with a canonical kinetic term in the Lagrangian, that describe dark matter 
halos in terms of bosonic scalar flelds, see e.g. Refs. [128, 66, 67, 68, 21, 63, 64, 65, 129, 130]. 
In this Section we assume that our scalar fleld configurations only depend on the radial 
direction. Three main results are achieved. First, we are able to find a purely kinetic 
Lagrangian which allows simultaneously to provide a flat rotation curve and to realize a 
unifled model of dark matter and dark energy on cosmological scales. Second, an invari- 
ance property of the expression for the halo rotation curve is found. This allows to obtain 
purely kinetic Lagrangians that reproduce the same rotation curves that are obtained start- 
ing from a given density profile within the standard Cold Dark Matter (CDM) paradigm. 
Finally, we consider a more general class of models with non-purely kinetic Lagrangians. 
In this case one can extend to the static and spherically symmetric spacetime metric the 
procedure used in Ref. [33] to find UDM solutions in a cosmological setting. Such a proce- 
dure requires that the Lagrangain is constant along the classical trajectories; we are thus 
able to provide the conditions to obtain reasonable rotation curves within a UDM model 
of the type discussed in Ref. [33]. 

9.1 Static solutions in Unified Dark Matter models 

Let us consider a scalar fleld which is static and spatially inhomogeneous, i.e. such that 
X < 0. In this situation the energy-momentum tensor is not described by a perfect fluid 
and its stress energy-momentum tensor reads 

= {P\\ + p)nt,nu - pg^u (9.1) 

where 

p = -p^ = -C, (9.2) 

riij, = V^yj/V— 2X and = C—2XdC/dX. In particular, is the pressure in the direction 
parallel to n^, whereas p± is the pressure in the direction orthogonal to n^. It is simpler to 
work with a new deflnition of X. Indeed, deflning X = — % we have 

= V^<^/(2x)^/' (9.3) 
n=2x^-p. (9.4) 

Let us consider for simplicity the general static spherically symmetric spacetime metric i.e. 

ds'^ = - cxp (2a(r)) dt'^ + exp (2/3(r)) dr^ + r'^dn^ , (9.5) 

where dil^ = d9'^ + sin^ 9d(j)'^ and a and /3 are two functions that only depend upon r. 
As the authors of Refs. [123, 95] have shown, it is easy to see that the non-diagonal term 
T*"* vanishes. Therefore if could be either strictly static or depend only on time. In this 
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section we study the solutions where depends on the radius only. 



In the following we will consider some cases where the baryonic content is not negligible 
in the halo. In this case we will assume that most of the baryons are concentrated within 
a radius rj,. If we define M* as the entire mass of the baryonic component then for r > ri, 
we can simply assume that M* is concentrated in the center of the halo. 

Considering, therefore, the halo for r > rf,, starting from the Einstein's equations and 
the covariant conservation of the stress-energy (or from the equation of motion of the scalar 
field, Eq. (2.4)), we obtain 

[r exp (-2^)]'} =p^^= Anpr^ , (9.6) 



dr 

-2 (exp [-2(a + P)] [r exp (2a)]' - l} = p,, ^ a' = J- ^J^^ ^ , 

(9.7) 



[■'■ A%r 

5iPt±i±M{[,expa]'/J'- [r(expa)']'} =P, (9.8) 



|f = -(P|,+p) M 

(which are the 00, rr and 69 components of Einstein's equations and the r component of 
the continuity equation respectively) where 

exp(-2/3(r)) = 1 - (M + M,)/(47rr) 

and 

R = ln[r^ exp(a(r))] 
, where a prime indicates differentiation with respect to the radius r. 



A first comment is in order here. If /3' = and ii) [r (exp a)']' > 0, then we can 
immediately see that p < 0. These conditions must therefore be avoided when trying to 
find a reasonable rotation curve. For example, neglecting the baryonic mass, the special 
case of p = A/r'^ and exp(Q;) ~ r"*, where A and m are constants, fall into this case. One 
thus recovers the no-go theorem derived in Ref. [95] under the assumption that the rotation 
curve ^ 1 is constant for all r. 

The value of the circular velocity Vc is determined by the assumption that a massive test 
particle is also located at 6 = 7r/2. We define as massive test particle the object that sends 
out a luminous signal to the observer who is considered to be stationary and far away from 
the halo. 

Considering the motion of a massive test particle, say a star, in a such a halo, its 
trajectory is then described by a curve = (t, r, 6, (p) parameterized by some affinc pa- 
rameter; here wc use its proper time r. Its four velocity is then simply = dx^/dr. Due 
to spherical symmetry, we can assume without loss of generality that the star's ecliptic is 
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located in the 6 = tt/2 plane. Since the star is a massive particle, its norm is u^/u^ = — 1, 
which becomes the constraint equation 

exp {2a)P - exp (2^)r2 _ ^2^2 ^ ^ ^g^^^Q^ 

where a dot denotes a derivative with respect to proper time r. Since the metric does not 
explicitly depend on 9, the star's angular momentum I is conserved, 

l = r'^4). (9.11) 

Similarly, the metric does not explicitly depend on t, and there is a conserved energy E, 

E = e^p{2a)i. (9.12) 

Substituting equations (9.12) and (9.11) into equation (9.10), one finds a first integral of 
motion for the star, 

If2 + V(r) = 0, (9.13) 

where its effective potential is 

V(r) = I exp (-2/3) (l + l--]- ^E^ exp [-2(a + /?)]. (9.14) 



Note that the potential explicitly depends on the energy. Stationary orbits at radius r 
exist if V and dV/dr vanish at that radius. The former condition yields 



;2 

1 + ^ = £;^exp[-2a(r)], (9.15) 



whereas the latter gives us 



/ 72 \ ;2 

-P' I 1 + - ^ + E^a' + 13') exp [-2a(r)] = 0. (9.16) 

Substituting equation (9.15) into equation (9.16) and using the Eqs. (9.6), (9.7) and (2.4), 
we get the following equation 

f/T^ _ (M + M,)(r)^^2^^ 



1 + Z7r2 87rr 

which directly relates the angular momentum / to the density profile of the halo. 

In this case, through the definition of the star's angular momentum I and Eq. (9.17), 
the value oi Vc = l/r can be rewritten as 

^2 ^ P||rV2 + (M + M.)/(87rr) 
" 1- [;j||r2/2 + (M + M,)/(87rr)] ' ^' ' 

but when we consider the weak-field limit condition (M + M*)/(87rr) <C 1 and since the 
rotation velocities of the halo of a spiral galaxy are typically non-relativistic, <C 1, 
Eq. (9.18) simplifies to [123] 

n M + M* piir^ 
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A second comment follows from the fact that the pressure is not small compared to 
the mass-energy density. In other words we do not require that general relativity reduces 
to Newtonian gravity (see also Refs. [124, 125, 126, 127]). Notice also that in the region 
where Vc ~ const. ^ 1 it is easy to see that in general exp(a) const, since from Eqs. (9.7) 
and (9.19) one obtains ro! k. v^. 

Finally, let us point out one of the main results (see also Ref. [38]). We can see that 
the relation (9.19) is invariant under the following transformation 

P — ^p = p + u{r) p\\ — >p\\=p\\+q{r) (9.20) 

if 

3g(r) + rq{r)' = -a{r) , (9.21) 

up to a proper choice of some integration constants. Thanks to this transformation we 
can consider an ensemble of solutions that have the same rotation curve. Obviously, these 
solutions have to satisfy Einstein's equations (9.6), (9.7) and (9.8), and the covariant con- 
servation of the stress-energy (9.9). Moreover, we will require the validity of the weak 
energy conditions, p > and p\\+ p'> 0, i.e. 

2exp(-2^ (9.22) 

r ox 

9.2 Unified Dark Matter models with purely kinetic Lagrangians 

Let us consider a scalar field Lagrangian C with a non-canonical kinetic term that depends 
only on X or x- Moreover, in this section we assume that = (or M ^ M*). 

First of all we must impose that C is negative when X < 0, so that the energy density 
is positive. Therefore, we define a new positive function 

gsix) = -C{X) . (9.23) 

As shown in Ref. [123], when the equation of state p|| = P\\{p) is known, one can write 
the purely kinetic Lagrangian that describes this dark fluid with the help of Eqs. (9.2) 
and (9.4). Alternatively, using (9.9), one can connect p|| and p in terms of r through the 
variable R. Moreover, it is easy to see that starting from the field equation of motion (2.4), 
there exists another relation that connects % (i-S- X) with r. This relation is 



X 



dx . 



^ (9.24) 



[r^ exp a{r)]^ 

with k a positive constant. If we add an additive constant to gs{x)i the solution (9.24) 
remains unchanged. One can see this also through Eq. (9.9). Indeed, using Eqs. (9.2) and 
(9.4) one immediately finds that Eq. (9.9) is invariant under the transformation p ^ p + K 
P\\ ~^ P\\ ~ K. In this way we can add the cosmological constant K = kio the Lagrangian 
and we can describe the dark matter and the cosmological constant-like dark energy as a 
single dark fluid i.e. as Unified Dark Matter (UDM). 

Let us notice that one can adopt two approaches to find reasonable rotation curves 
Vc{r). A static solution can be studied in two possible ways: 
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i) The first approach consists simply in adopting directly a Langrangian that provides a 
viable cosmological UDM model and exploring what are the conditions under which 
it can give a static solution with a rotation curve that is fiat at large radii. This 
prescription has been already applied, for example, in Ref. [123]. 

ii) A second approach consists in exploiting the invariance property of Eq. (9.19), with 
respect to the transformation (9.20) (when the condition (9.21) is satisfied). Usually 
in the literature one reduces the problem to the Newtonian gravity limit, because 
one makes use of a CDM density profile, i.e. one assumes that in Eq. (9.19), py ^ 
M/(47rr^). We can therefore use Eqs. (9.20) and (9.21) to obtain energy density 
and pressure profiles p{r) and p\\{r) that reproduce the same rotation curve in a 
model with non-negligible pressure. Next, we find an acceptable equation of state 
P\\ = P\\{p) such that we can reconstruct, through Eqs. (9.2) and (9.4), the expression 
for the Lagrangian C Such a procedure establishes a mapping between UDM and 
CDM solutions that predict the same halo rotation curve Vc{r). As a starting point 
we could, of course, use very different CDM density profiles to this aim, such as 
the modified isothermal-law profile [131], the Burkert profile [132], the Moore profile 
[133], the Navarro- Frenk- White profile [134, 135] or the profile proposed by Salucci 
et al. (see for example [136, 137, 138, 139]). 

As we have already mentioned, the possible solutions one finds in this way have to 
satisfy the Einstein equations (9.6), (9.7) and (9.8), the conservation of stress-energy 
(9.9) and the weak energy conditions. Moreover, the resulting UDM scalar field 
Lagrangian must be able to provide cosmological solutions that yield an acceptable 
description of the cosmological background (sec, e.g., Ref. [33]) and low effective 
speed of sound (see for example Refs. [70, 100, 32]) so that cosmic structure formation 
successfully takes place and CMB anisotropies fit the observed pattern [27, 31, 25, 28). 

Below, using approach i), we provide a worked example of a UDM model with purely 
kinetic Lagrangian which is able to describe a fiat halo rotation curve and then, using 
approach ii), we give a general systematic procedure to obtain a possible Lagrangian of 
UDM model starting from a given CDM density profile. 

9.2.1 Approach i): The generalized Scherrer solution 

Let us consider the generalized Scherrer solution models obtained in Ref. [33] (see also 
Sect. 4.3). These models are described by the following Lagrangian 

C = -k + gn(x-xY (9.25) 

where > is a suitable constant and n > 1. The case n = 2 corresponds to the 
unified model proposed by Scherrer [30]. If we impose that today [{X — X)/X]" ^ 1, the 
background energy density can be written as 

p{a{t)) =pK + PDM , (9.26) 
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where p\ behaves hke a "dark energy" component (pA = const.) and pDM behaves hke a 
"dark matter" component i.e. pDM oc a~^, with a{t) the scale factor. 

A static solution for the generalized Scherrer model can be obtained in two possible ways: 

1) Starting from the analysis of Ref. [95], in the case of a barotropic Lagrangian for 
the homogeneous field. The authors of Ref. [95] indeed concluded that for n ^ 1 
flat halo rotation curves can be obtained. In particular they studied spherically 
symmetric solutions with the following metric, 

ds^ = -(^y^ dt'^ +N{r)dr^ + r^dn\ (9.27) 

where r^^ is a suitable length-scale and b = 2v'^. In the trivial case where N{r) is 
constant they find JC{X) oc X"^/^ with 6^1. For X » X the Lagrangian £ = 
—A + gn{X — X)^ takes precisely this form. 

2) In the analysis of Ref. [123], solutions where cp is only a function of the radius are 
considered. When the Lagrangian has the form £ oc X", with n ~ 10^ the halo 
rotation curve becomes flat at large radii. In this case n must be an odd natural 
number, such that the energy density is positive. Our model is able to reproduce this 
situation when the matter density is large, i.e. when |X| S> X. 

Alternatively, if we wish to avoid large n (c.f. case 2) above) we can start from the following 
Lagrangian 

C = -A + exgn{\X\-Xy (9.28) 

where ex is some differentiable function of X that is 1 when X > X and —1 when X < 
—X < 0. In this way when X > X > we recover the Lagrangian of the generalized 
Scherrer solutions. When X < and x = > X we get 

£. = -A-gn{x-Xy (9.29) 

and, with the help of Eqs. (9.2) and (9.4), we obtain 

p = -p± = -C, p\\ = (2n - l)gn (x - ^) " + 2ngnX (x - x) - A . (9.30) 

Now, requiring that x be close to X (i.e. ^x ~ ^) 2ngnX ^x ~ S> 0{A), 

and starting from the relation (9.24) that connects x with r, we get 

^ r2 exp (a(r)) " ^^"^^^ 

Consistency with our approximations implies that we have to consider the following ex- 
pressions for radial configurations with r bigger than a minimum radius r^m- Iii this case 
P|| and p become 

^11 ~ r2exp {air)) ' ^ ~ [r2 exp (a(r))]V(n-i) (^"^2) 
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where A = 2{kXfl^ and B = g^, k^l"^ jing^X^I'^) 

Using Eqs. (9.6) and (9.7), one can calculate the values of the metric terms exp (a) and 
exp (/?) and, thus the value of p and Alternatively we know that when Vc ~ const. <C 1 
at large radii, in a first approximation, we can set exp(a(r)) « C = const. Therefore for 
n 7^ 3, we can write the function M as 

47r B / n — 1 n-a \ 
^(0 ~ —77^ 1?— + D (9.33) 

where we could also set D = for n > 3. Instead, when 1 < n < 3, the second term has to 
be larger than the first one. 
In these cases becomes 

2. - A , B fn-l 1 D 



For n = 3 we have 



^^^^^ " 2C + 2 CV(n-l) L-3rV(n-l) + 7 ' • ^^-^^^ 



47r B / r \ 



where r > r and 



9 -4 S 1 /r\ M(f) 

In other words we see that the circular velocity becomes approximately constant for suffi- 
ciently large r. 

However, let us stress that exp(a(r)) cannot be strictly constant, and that it should 
be chosen in such a way that the positivity of Eq. (9.8) is ensured. 

This example can be generalized also to M* ^ 0. Obviously, in such a case we have to 
assume that r > r;, > rmin- In this case A;, rmin-, A, B (through exp(/3(r))) and C depend 
on M*. 

The spherical top-hat solution for this model, which provides the link with the cosmo- 
logical initial conditions, is described in the Sect. A. 

9.3 Approach ii): A general prescription to obtain UDM Lagrangians starting 
from a profile of an energy density distribution of CDM 

Defining the energy density distribution of CDM as pcdm(?') (with pcDM = 0), the trans- 
formation (9.20) becomes 

p(r-) = pcDuir) + C7{r) , (r) = q{r) . (9.37) 

Now, starting from a given CDM density profile, through Eqs. (9.6), (9.7), (9.9) and (9.21) 
we can determine exp (a) , exp(/3), p and py. In a second step we provide the conditions 
to ensure that the energy density is positive ^. In this case, after some simple but lengthy 



^Thanks to this condition, through Einstein's Eq. (9.8), we can evade the no-go theorem derived in 
Ref. [95]. 
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calculations, one finds 



2'W {r^^^-2rQir)^ - 
+Q(r) f4r + 3%^+4rVDM 



2Q\r) 

rMcDM(r) 



47r 



(4 + SrVcDM) 



i3(r) 
Air) 

a{r) 



Q{r) 

Q(r) + g(r) 

2H(r) 
1 - Q'{r) 



McDuir) 



4tt 



(9.38) 

(9.39) 

(9.40) 

(9.41) 



where Q(r) = r{r'^q + 1) , B(r) = rexp(— 2/3) and A{r) = (ra' + 1). Here we define 
Mcdm('') = 47r [q r'^pcDuir) dr. At this point it is easy to see that Eq. (9.38) does not 
admit a simple analytical solution for a generic pcDM- On the other hand we know that, 
through pcDM, all these functions depend on the velocity rotation curve Vc{r). Moreover 
v'^{r) ^ 1. Therefore, defining Vc as the value that Vc assumes when the rotation curve is 
fiat at large radii or the maximum value of Vc with a particular profile of pcDM) one can 
expand Q, A and B as 



Qir) = Q(0)(r) + t^c'Q(i)(0 + ^Q(2)(r) + • • • , 
Air) = ^(o)(r) + ^;^^(i)(r) + %^^(2)(r) + . . . , 



2! 

Bir) = i3(o)(r) + v',B^,)ir) + -^%)(r) + . . . . 



(9.42) 



Following this procedure one can determine p and in a perturbation way, i.e. 



2^2 



(^c) 

Pir) = p(o){r) + v'^p^i){r) + ^|-p(2)(r) + . . . , 

(^2)2 

P\\ (r) = p\\ (0) (r) + vIp\\ (1) (r) + -^P\\ (2) (r) + ■ ■ 



(9.43) 



(9.44) 



Now, looking at the various CDM density profiles which have been proposed in the literature 
[131, 132, 133, 134, 135, 139], we see that one can always take pcDM as 

pCBuir) = VcPcBM ii){r) , 

then 

McBuir) = v^McDM (i)ir) = in / fVcDM (i)(^) df . 

Jo 

For the zeroth-order terms one immediately obtains 



(9.45) 
(9.46) 



2(0) 


= r 


-4(0) 


= 1 


^(0) 


= r 



(9.47) 
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At the first order one gets 



2 

2(1) = - / r^PCBM{i){r) dr , 
f Jo 

1 AfcDM (!)(?') 



A, 



(1) 



2r Att 

- r pcDM (1) (r) . (9.48) 



For completeness we write also the second order for Q 

1 r ,,^CDM(l)(r) 

r Jo 



1 r ,_MCDM(1)( 
2(2) = - / 

^ ' r Jo 4tt 



Q(l)(f) -fVcDM (!)(?') 



(9.49) 



Let us stress that if one considers also terms O {vf), Eq. (9.18) instead of Eq. (9.19) should 
be used. In such a case, Vc slightly changes with respect to the velocity rotation curve that 
one obtains using a CDM density profile. 

For our purposes we can consider only the zeroth and the first-order terms. At this 
point, one can finally calculate the value of p and One gets 

p{r) = pGHMir) + = ^l\-A r^PCDM{i){r) df - PcDM(i){r)j , (9.50) 

As far as the values of the metric terms exp (a) and exp [j3) are concerned, we obtain the 
following expressions 



exp {2a) = exp (2a(f)) exp 



_2 / l^CDM(i)(r) 
Air 



(9.52) 



exp (-2/3) = 1 + ^ (^2 fVcDM (1) (r) df - ^ ^CDMJi)(r) ^ _ 

Now, it is immediate to see that if we want a positive energy density we have to im- 
pose 2 /J r^PcDM (1) (^) df > Pcdm (!)('')• Prom Eq. (9.6) we know that M(r) = 
47r /r^^ f^p{f) + M(r(o)) and McDM(r) = Att J^f'^pcDM{f) df + Mcdm(^). Therefore we 
need to know what is the relation between f and r(o)- This condition is easily obtained if 
we make use of Eq. (9.19). Indeed, we get 

M(i)(f(o))-McDM(i)(r) , 2 r(o).3 r(o),2 

—^-^-^ — '-^ ^ — / ^ PCDM(i){r) dr= r PcDM(i)(n dr , (9.54) 

471" ^(0) Jo Jf 

which finally guarantees the invariance of the rotation velocity with respect to the trans- 
formation in Eqs. (9.20) and (9.21). 

Let us, to a first approximation, parametrize the various CDM density profiles, at very 
large radii (i.e. when we can completely neglect the baryonic component) as 

PCDM = (9.55) 
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where k is a proper positive constant which depends on the particular profile that is chosen 
[131, 132, 133, 134, 135, 139]. For example for many of the density profiles the slope is 
n = 3 for large radh [132, 133, 134, 135, 139]. 

In this case a positive energy density p > requires n > 2. At this point let us focus 
on the case where 2 < n < 4, since this gives rise to the typical slope of most of the density 
profiles studied in the literature. Therefore we obtain for p{r) and p\\{r): 



f \ —0 Th ^ \ . , _o 2 1 , . 

p(r] = v„ n , Oil (r) = v^. k- . (9.56) 

In particular, 

1) for n = 2, we get 

p{r) = 0, p||(r) = PCDM = Vc , (9.57) 

and for the relation between f(o) and f one can choose, for example, f(o) = f = 0. In 
other words, for large radii we have that p(r) <C P\\{r)- 

2) Also for 2 < n < 3 one can choose f(o) = f = 0. 



r.2 1 



3) For = 3 



2 

p(r) = PCDM, P\\ir) =vl K-^ , (9.58) 



and, through Eq. (9.54), we have to impose that 



M(i)(f(o)) - McDM (i)(r-) _ 1^ ^M', _ 2 . (9.59) 



47r \ r 

Notice that the energy density profile is the same as the CDM one only for large radii 
so that M(i)(r) differs from Mcdm {Vjir)- 

4) In addition, for 3 < n < 4, also through Eq. (9.54), we have to impose that 

M(i)(f(o)) - Mcdm (i)(r) _ f^-" (n - 2) 3_„ 

47r n-3 (4-n)(n-3) W " ^ ' 

Now let us focus where 2 < n < 4. Starting from Eq. (9.56) to express py = p\\ (p) we solve 
Eq. (9.4) to recover the Lagrangian for the scalar field 

n 2k n 

p(x) = -C = kx^^ , p{x) = 7 ^X^^ (9.61) 

(n-2) 

where is a positive integration constant. We can see that, for this range of n, the exponent 
is larger than 1; thus there are no problems with a possible instability of the Lagrangian 
(see Refs. [123, 78, 77]). Therefore, through the transformation p — > p + A p|| — )• p|| — A, 
this Lagrangian can be extended to describe a unified model of dark matter and dark 
energy. Indeed, starting from the Lagrangian of the type (9.28), when |X| » X and if 
k = Qn, C takes precisely the form (9.61). 

Finally, we want to stress that this prescription does not apply only to the case of an 
adiabatic fluid, such as the one provided by scalar held with a purely kinetic Lagrangian, 
but it can be also used for more general Lagrangians C{ip,X). 
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10. Conclusions 



In this work we explored the possibihty that the dynamics of a single scalar field can account 
for a unified description of the Dark Matter and Dark Energy sectors, leading to a Unified 
Dark Matter (UDM) model. In comparison with the standard DM + DE models (e.g. 
even the simplest model, with DM and a cosmological constant), in UDM models there are 
two simple but important aspects to consider: first, the fluid which triggers the accelerated 
expansion at late times is also the one which has to cluster in order to produce the structures 
we see today. Second, from the last scattering to the present epoch, the energy density of 
the Universe is dominated by a single dark fluid, and therefore the gravitational potential 
evolution is determined by the background and perturbation evolution of just such a fluid. 
As a result the general trend is that the possible appearance of a sound speed signiflcantly 
different from zero at late times corresponds to the appearance of a Jeans length (or a sound 
horizon) under which the dark fluid does not cluster any more, causing a strong evolution 
in time of the gravitational potential (which starts to oscillate and decay). Specifically in 
this paper we have explored UDM models defined by Lagrangian of k- essence models. This 
allows to find suitable solutions around which the scalar field describes a mixture of Dark 
Matter and Dark Energy. Finally we also investigated the static and spherically symmetric 
solutions of Einstein's equations for a scalar field with non-canonical kinetic term. 
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A. Spherical collapse for generalized Scherrer solution models 

Let we assume a flat, homogeneous Priedmann-Robertson- Walker background metric. In 
such a case, the background evolution of the Universe is characterized completely by the 
following equations 

H = -\{P + P). (A.2) 
where the dot denotes differentiation w.r.t. the cosmic time t. 

Now let us consider a top-hat spherical over-density with the purely kinetic model with 
the Lagrangian C = —A + gn{X — X)" and with Qn > 0. For this particular case within 
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the over-dense region we have a single dark fluid undergoing spherical collapse, which is 
described by the following equation 



{PR + 3pii) 



(A.3) 



R _ 1 

where R, pR and pR are respectively the scale-factor, pressure and energy density of the 
over-dense region; pR and pn are defined by the following expressions 

PR = K + 2ngnX{XR - X^-^ + (2n - l)gn{XR - X)^ 



PR = 9R = -^ + 9c{Xr - X)" 

with Xr = X{R) a function of time. 
The equation of motion is 

(dg_ 



dXR 



dXf,J dNR 

where (INr = dR/R. The solution of Eq. (A.6) (for dgR/dXR,XR ^ ) is 



Xr 



(A.4) 
(A.5) 



(A.6) 



(A.7) 



where we can choose kR = 
Replacing Eq. (A.5) in Eq. (A.7) 

Xr 



we find 



, with Rta the value of R at turnaround. 



- ta 



ngn{XR - X) 



n—l 



kRR- 



(A.8) 



Using now the explicit expressions for pR and pR we arrive at the following set of equations 

R _ 1 
^ ~ ~3 

tR~ 



ri^gl 



^A + ngnX{XR - X^'^ + (n + l)gn{XR - X)" 
{Xr - Xf^-' + X{Xr - X)2("-i) 

For {Xr - i')/i: < 1 Eq. (A.9) becomes 

-A + ngn\XR,^ - Xr\XR,^xf2 



(A.9) 
(A.IO) 



R 
R 



R 

Rta 



-3 



We can now write all the equations that describe the spherical collapse 

2 



a/ 6 
PA = A 



{pA + Pdm) 



PDM = 2ng„\Xt^ - X\^-'{XtaX)-2 



a 

ata 



R 
R 



= --^ {PRdm - '^PRa 



PR^^ = 2ngr,\Xta-Xr-\XR,^X)-2 



R 

Rta 



(A.ll) 

(A.12) 

(A.13) 
(A.14) 

(A.15) 
(A.16) 
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where ata = a{tta)- 

Following now the same procedure of Ref. [140] we can define x and y 

x = — (A.17) 

ata 

y = A . (A.18) 
In this way we can redefine /9dm and PRdm such that 

X"^ 

PRdm - C ^3 (A.2U) 

where 0dm is the (k-essence) dark matter density parameter, and ^ = (p/pDM)|a;=i- Then 
Eqs. (A. 12) and (A. 15) become 

fiT 1 

^ = (xSJ„„(i))-' , (A.21) 

= -^[<-2,»if.]. (A.22) 



where dr = Hta^/n^^j{x = 1) and Ka = /?^/[3i7i2^Jl^^,(z = 1)]. 

Defining U as the potential energy of the over-density and using energy conservation 
between virialization and turnaround, 



'rr 

^+2aS 



Uta , (A.24) 



-I vir 



we obtain 

(1 + q)y - 2qy^ = ^ (A.25) 



where 



, = ('^) = ^ , (A.26) 

Pjy=l C 



in full agreement with Ref. [141]. 
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